CHEMICAL ENGINEERING SCIENCE 
GENIE CHIMIQUE 


OCTOBER 1955 








Gas mixing in long pipelines 


J. F. Davinson, D. C. Farqunarson, J. Q. Picken and D. C. Taytor 
Department of Chemical Engineering, Tennis Court Road, Cambridge, England 


( Received 25 March 1955) 


Summary—tThe longitudinal diffusion coefficient in a 8 in. diameter pipe carrying refinery 
gas over a 9,000 ft distance has been measured by using hydrogen as a tracer. The results 
show that the formula for the effective diffusion coefficient given by Taytor, K = 7-14 aU \/cy, 
may be used for gases as well as for liquids. This means that the Schmidt group has no influence 
on longitudinal diffusion. 


Résumé—En utilisant de lhydrogéne comme indicateur, les auteurs ont mesurés le coefficient 
de diffusion longitudinal dans une conduite de 3° (7,62 cm) servant au transport des gaz de 
raffinerie, sur 9.000 ft (2,743 km). Les résultats prouvent que la formule de TayLor pour le 
coefficient de diffusion effective, K = 7,14 aU cf est valable pour les gaz comme pour les 
liquides. En conséquence, le groupe de Scumrprt est sans influence sur la diffusion longitudinale. 


Wuewn fluid flows through a pipe with turbulent 
motion, the dispersion of matter can be predicted 
from a simple formula given by Tay.or [5]. He 
has shown that the longitudinal transfer, relative 
to a plane moving with the mean velocity of flow 
in the pipe, follows the ordinary diffusion equation, 
and the virtual coefficient of diffusion is 
K =714aU yc. a is the radius of the pipe, U 
the mean flow velocity and c¢, the friction co- 
efficient. 

The object of the present work was to find out 
whether this result would be influenced by a wide 
variation in the Schmidt group y/pD. Such a 
variation can be achieved by changing from 
liquids, which were used in all the experiments 
quoted by Tay.or, to gases. Experiments on 
mass transfer from a solid boundary to a turbulent 
stream have shown [4] that the Schmidt group 
has an appreciable effect, and a similar result 
might be expected with longitudinal dispersion 
in pipes. The experiments described here deal 
with a 9,000 ft stretch of 3 in. diameter pipeline 
carrying refinery gas, and the Schmidt group was 
about 1,000 times smaller than that of previous 
investigators. 


The theory given by Taytor refers to a 
straight pipe, but in the present case the pipe 
contained a number of sharp bends. Calculations 
showed that the bends would contribute only 5% 
of the pressure drop, and it was therefore antici- 
pated that the effect upon the mixing would be 
small. Part of the interest of the experiments lay 
in finding out whether this expectation would be 


justified. 
GENERAL PROCEDURE 


Following the method of previous workers, a small 
shot of tracer material was injected at the 
beginning of the pipeline. The material then 
moved down the line with the mean velocity, 
about 8 m.p.h., spreading longitudinally and 
giving a measurable concentration at the outlet 
for about 30 seconds. The tracer was hydrogen, 
and a conductivity meter was used for the 
analysis. The lag of the instrument was such that 
it would detect the arrival of the injected slug, 
but would not give an accurate concentration- 
time graph. Two shots of hydrogen were therefore 
injected, with a 4 minute interval. The arrival 
of the first shot, about 15 minutes after injection, 
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was detected by the meter connected to the Sampling 
sampling point. 8 minutes 45 seconds later the The sampling arrangement, shown in Fig. 1, was 
first of 12 evacuated sample bottles was connected — similar to that of Danckwenrts et al. [1]. Eva- 
to the line at the sampling point, and by taking cuated sample bottles of about 100 cc. capacity 
the other samples at 3 second intervals, the were connected to the line via a two way cock 
distribution of hydrogen near the time of arrival and a manifold. The sample bottles were at most 
of the second shot could be obtained by analysing — 6 ft from the pipeline, and by taking a continuous 
the samples at leisure. bleed of gas through the manifold and the two 
The pipeline worked at about 200 p.s.i. gauge way cock at about 20 ft ‘sec, a lag of only about 
with a pressure drop of about 8 p.s.i. over the 4 sec between the line and sample bottle com- 
test length, and the latter quantity had to be positions was obtained. 
known accurately for the calculation of c,. The 


. ming 
method used was to connect simultaneously a Timing 


cylinder to each end of the line, and subsequently On turning the two way cock, a thread of mercury 


to measure the pressure difference between the rose in the evacuated manometer, and by moment- 


cylinders on a differential pressure gauge. arily joining the platinum wires, made a mark on 
a uniformly rotating drum, At first it was found 
that the capillary tube fractured at A (Fig. 1), 


EXPERIMENTAL DETAILS Py. 
owing to the high velocity of the mercury. This 


Hydrogen injection difficulty was overcome by putting a small piece 
Following the method employed by DancKWerTs of asbestos wool inside the tube at A. 
et al. [1], a cylinder, fitted with a pressure gauge The marks on the drum fixed the time intervals 
and connected to the line via a quick acting cock, between the samples, and the time of opening 
was filled with hydrogen from a storage cylinder. the first tap was noted, from a watch, to the 
From the volume of the intermediate vessel nearest second. In this way the sampling was 
(0-170 cu. ft.), the pressure drop on injection co-ordinated with the injection, which had been 
(from 400 p.s.i.g. to 260 p.s.i.g.), and the tem- made at an agreed time. 
perature, it was possible to calculate the quantity 
of hydrogen injected. Analysis 
The Cambridge differential katharometer [2] 
already employed to detect the arrival of the 
signal pulse, was adapted to measure the con- 
centration of hydrogen in the sample bottles. 
When the instrument is operating normally, heat 
is conducted through the gas from an electrically 
heated wire to the walls of a small chamber in 
a metal block. The gas is supplied to the meas- 
uring chamber by diffusion through small drillings 
connected to a hole carrying a continuous flow 
of gas. The instrument is fitted with two similar 
LEED <> chambers for the comparison of gases. The 
heated wire in each chamber forms an arm of a 
WI Wheatstone bridge, and the heating current also 
serves as a bridge input. 
For the analysis of samples, the block was 
fitted with glass taps so that the measuring 
ecaaeoeed chamber could be evacuated, and subsequently 
Fic. 1. Sampling arrangement. connected to the sample bottle. It was found 
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that the instrument needed many minutes to give 
a steady reading after admitting a sample. The 
time was reduced to about 8 minutes by con- 
necting fresh bottles to both sides of the instru- 
ment at once. One bottle contained the sample, 
and the other contained gas from a cylinder which 
had been filled from the pipeline just before 
taking the samples. 

The instrument was calibrated by mixing a 
known volume (about 100 cc) of line gas with a 
measured volume (about 2 cc) of hydrogen, giving 
readings of the same order as those in the tests. 

A daily gas analysis, using a mass-spectrometer, 
was done by the refinery, and from these results 
the density was calculated. On a molecular 
basis, the gas mainly 
methane, ethane and ethylene, as well as small 
amounts of propane, propylene and _ heavier 
hydrocarbons. There were daily variations in 
hydrogen content of about 10%. This variation 
affected the calibration of the katharometer, but 
did not influence the calculation of the diffusion co- 
efficient K, which depended only on the relative 
readings of the instrument on any single run. 


contained hydrogen, 


Measurement of pressure drop 


This was done by means of a Budenberg “Micro- 
var” differential pressure gauge reading differ- 
ences up to 80 lb./sq. in. The cylinders A and E 
Fig. 2, were filled at opposite ends of the pipeline. 
With B closed, C and D were connected together, 
and the pressure difference was read from the 
gauge G, with a small correction for the volume 


of pipework between B and C. 


pA 
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Fic. 2. Pressure drop measurement. 


RESULTS 
Friction coefficient 
Table 1 gives values of the friction coefficient Ch 
calculated from the pressure drop and from the 
simultaneous measurement of the mean velocity 
by the injection and sampling technique. 


Table 1 





Friction 
coefficient 


108 of 


Pressure 
drop 
p.8.i. 


Density | 
lb./eu. ft. | 
0-73 . 2. 5-2 
077m 8 «(| 2. 5-2 
0-74 . 3-4 54 





The values are roughly in accordance with the 
graph given by GoLpsTern [3] for a roughness to 
pipe diameter ratio of 0-0015. This graph was 
used for estimating the friction coefficient at other 


values of R. 


Parameter K /aU v/cy 

The katharometer readings, giving a measure of 
the hydrogen are 
plotted in Fig. 3. Also shown is the function 
c = Aexp — B(t — t,)*, the values of A, B and 
t, being determined by the method of least 


change in concentration, 


squares. An alternative method is to treat t as 
the variate, and ¢ as the frequency, so that 
to = Le, t,/ Ze, 1/2B = Let?/ Ze, — ty’. 
This method is less accurate, but its advantage is 
that the standard deviation of B may be calculated 
without undue labour. 

On applying these methods to calculate B, and 
getting the effective diffusion, coefficient K from 
the relation K = L?,/+ *B, we obtain the results 
shown in Table 2. L is the length of the pipe, 
9,000 ft. 


and 


Table 2. 





K/aU Vc, from: 
Least Fre- 
squares 


























Volume injected 
From the areas under the curves shown in Fig. 3, 


and the katharometer calibration, the volume of 
hydrogen injected in each case can be estimated. 





KATHAROMETER DEFLEXION —e= 
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The results were 8-10°%, greater than the values 
given by the pressure drop at the injection point. 
This discrepancy may be due to variation in the 
pipe diameter which could not, of course, be 
measured. 
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Fic. 3. Concentration of hydrogen at the outlet. The 
curves are fitted by the method of least squares. 


DiscuSssION 
The theoretical value of K (aU yc, = 7-14, given 
by Taytor [5] and previously verified only for 
liquids, is appreciably lower than the values 
shown in Table 2. However, the difference is of 
the same order as the standard deviations and we 
cannot conclude that there is a definite difference 
between gases and liquids, in view of the following 


uncertainties : 


(i) The effect of the bends, 20 elbows at 90° and 
8 at 45°, is uncertain. Although these bends 
account for only 5% of the pressure drop, it is 
impossible to be sure of the magnitude or the 
direction of their contribution to the longitudinal 
dispersion. 


(ii) It is possible that some initial mixing may 
have been caused by the method of injection. The 
hydrogen was injected with the velocity of sound 
through a side tube, and this might distribute the 
gas over (say) a 20 ft. length instead of the 2 ft. 
assumed theoretically. The effective length of 


the slug at outlet was roughly 100 ft., so that an 
initial spread of 20 ft. would not be negligible. 


CONCLUSIONS 


It is certain that the effect of Schmidt number on 
longitudinal mixing is very much smaller than its 
effect on transfer to the walls. In the present case, 
using the commonly accepted factor Sc’ [4], the 
wall transfer coefficient would be 100 times higher 
for gases than for liquids. The widest variation 
that could be inferred from our results would be 
a factor of 2, and it seems more likely that there 
is a negligible difference between gases and 
liquids. This conclusion seems reasonable, since 
transfer to the walls involves molecular diffusion 
through the relatively stagnant boundary layer, 
while longitudinal transfer is brought about by a 
combination of the gradient of the mean velocity 
with mixing due to eddies. 
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NOTATION 


peak deflexion of the katharometer 
radius of the pipe 
= decay constant in A exp — B(t — &)* 
concentration of hydrogen 
= pipe friction coefficient 
= molecular diffusivity 
effective longitudinal diffusivity 
= length of pipeline 
R = Reynolds number = 2pUa/pz 
Sec = Schmidt number »/pD 
t = time 
ys = viscosity of line gas 
p = density of line gas 
Subscripts 
0 = peak concentration 
i = number of experimental point 
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Summary 


In previous papers the computations, both theoretical and numerical, involved 


in rectification for multicomponent and complex mixtures in cascades, were considered. The 
formulae may be extended to performance in packed columns so that the calculation of the 
number of transfer units for a given separation may be obtained with ease. 


Résumé 


la rectification par plateaux des mélanges a plusieurs composants. 


Dans des articles précédents, les auteurs ont étudiés théoriquement et numériquement 


Les resultats peuvent étre 


étendus au cas des colonnes A garnissages de telle sorte que le calcul du nombre des unités de 


transfert pour une séparation donnée peut étre obtenu facilement. 


MULTICOMPONENT SyYSTEMS* 


Ir n is the number of transfer units and y,(t) 
and y,*(i) is the vapour composition and the 
vapour composition which would be in equili- 
brium with the liquid, respectively, then the 
differential equation which describes the perfor- 
mance of a packed distillation column is 


dy,(%) 


1 
dn Yai?) (1) 


Yn® (¢) 
As in previous parts the equilibrium relation 
will be assumed to be of the form 


a (t) 2,,(7) 


2 « (i) tali) 


Yat?) 


while the overall mass balance over the top n 
transfer units and total condenser is 


Tol) 


Rai (3) 


Yal?) - R i j Tal?) T 
UnpeErwoop [1] considered this system and his 
results will be rederived so that they may be 
more easily generalized to complex mixtures. 
Eqs. (1), (2) and (8) may be combined to give 


* Notation is as in Parts 1-4. See this journal. 


R dx,{i) 
R+1 dn 


a (i) 2,(7) 


2 « (i) 2,(i) 


R+1 


It would be desirable to have formulae for 
the number of transfer units as a function of 
composition or for the change in composition 
as a function of the number of transfer units. 
With these ends in view a solution to Eq. (4) 
must be found. 

Using the functions »,(#), 

pi) 

An — p(t) 

defined in Part 1, eq. (3-8), one may assume in 
order to find the solution that 2,(i) can be 
expanded into a finite sum. 


z,(i) = Toft) © Qu(n) va(i) (5) 


v,(i) = 


It follows that 


i . = 


roi) pi) J" 
P . 
i [Ae — oar 
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from the orthogonality property of the set of 
v,(i). See Part 1, eqs. (3-8)-(3-11). 
A summation of eq. (5) on ¢ produces : 


1 
y ; 
“ Qn) = 5 
so that eq. (4) may be written 
R yes 
R+1, dn 
a(t) 2 Qn) v, (7) 


k 
74,2 —=>... 
oe TS ah 


v,(t) 


yy To y %op(i) 
i Ax — p(t) i P(t) (Ax — pl?) Ak 
and 

y dQ,(n) . 


+dQ,(n) 
x adn “ 


dn leat) + 1) 


y (Se Ag Ue(t) 
k an pl?) 


(7) 


since 
sy 4Q),(n) 


zt adn 
Equating coefficients of v,(i) there results 
R dQ,(n) Q,(n) 


= un) ~ RQ An) (8) 
zt & 


Eq. (8) may be integrated between 0 and n to 
give 
n 
Qi(m) = 
Q.(0) Ay 


The integral is a function of n alone and may be 
written f(m) so that 


R |n 


A, In @(n) 

Q(0) 

Thus if one considers two different k’s say 
k and j, then the function f(n) need not be 
computed at all since one can easily show that 


+ nA = f(n) (9) 


Ak 


Aj 
1 n Q.( n | Qi 0) 


10 
Q)(n) ne) 


n= —___. ] 


Ai Ae 


&,(0) 


Note that this formula gives the number of 
transfer units in terms of readily calculable 
quantities. The Q,(n) are computed from eq. 
(7) and the A, from 


p(t) xof?) 
i Ay — p(t) 


R 4 
Q,(0) = 
Hence 
1 , 


In 
A; — »% T; 


r Mm gg 2@,(t) 
° B+iG¢k—0 


ComPpLeEX MIXTURES 
In complex mixtures it is clear that eqs. (1) 
and (3) still hold with 2,(i) and y,(i) as distribu- 
tion Eq. (2) requires appropriate 
modification with the summation replaced by an 
integral. The 
applied to the multicomponent case in the same 
fashion as was done in Part 1. It is clear that the 
functions p(A) and w(A) may be defined as before 


functions. 


limiting process may now be 


to give 
A Xo(A) 7 cot mp(A) = R — 


@ 

P| (2d) dq = w (A) 
A—q 

0 


(11) 
where 2,(A) is the distribution function of the 
overhead product. 

Eq. (6) and (7) then become on passing to the 
limit 
Q (A, n) = C(A) 


@ 


P| ve dg + a,(A) cot mp(A)] (12) 


0 


Ar9(A) 


at >) 4. fo) 2 (13) 
w*(A) + [A 2o(A)] 


C(a) = 


while from eq. (5) 
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z,(p) -_ prop) 


Note that Q(A, n) depends upon the distribution 
function at two points within the column, for 
z,(qg) occurs explicitly in eq. (12) while 2,(A) 
occurs in eq. (13). The limiting process used in 
Part 1 then may be extended to eq. (9) to give 


QA,n) 
Aln + nA = F(n) 
Q(A, 0) 
If two different values of A are chosen, say A 
and y, then 


1 n [rl [exe 4 
Q(A,0)} LQ(y.") 


Thus one may compute the number of transfer 
units required for a given separation provided the 
two compositions are known in detail. The problem 
would generally not be stated in this form and 
a trial and error procedure of the usual kind 
would be required for other cases. 

Eq. (14) may be simplified somewhat for 


n= (14) 


<—- 


Q(A,o) = C(A) 
p (gan dd + Q(A,n) = cot a) G 


P| ap dq + %9(A) » cot “A 


R +1 
= (2) = 


Therefore, if, 


QAn) _ 
Qo) © 


R+1 
to [S(A,n)}* 


y— A [S(yn)] 


S(A.n) = 


f “ioe dq + a,(A) m cot “) 


This formula is much better from the point of 
view of computation. It should be pointed out 
that the integrals appearing here although 
written from 0 to @ are not really improper at 
either end since relative volatilities are not 
zero or infinite. The integrals are improper 
at A = p, however, and the computations require 
the Cauchy principle values. The computation 
of such integrals was discussed in Part 8. 
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8.5. 


THE power consumption was measured as des- 
cribed in chapter 2.2. 

The abundance of data given in literature on 
the power consumption of stirrers made it super- 
fluous to investigate this point extensively. We, 
therefore, only measured the power consumption 
of stirrers under conditions not yet, or not 
adequately, dealt with in the literature. In this 
way various stirring conditions on a basis of 
mixing time and power consumption can be 
compared (chapter 4). 

The power consumption P was converted in 


the usual the group 
2 


age which was then correlated with the other 
p 


Power consumption 


way into dimensionless 


dimensionless groups, such as Reynolds number 


nd* , ; 
Pp - and various length ratios. These correlations 


be 
are of an experimental nature. It was not possible 


to establish a relationship between power con- 
sumption of the stirrer and mixing time. The 
two quantities are more or less independent of 
each other. Depending on stirring conditions, it 
is possible to obtain several mixing times at a 
given power consumption of the stirrer. Naturally, 
the shortest mixing time is generally to be 
preferred. 

Paddle stirrer—A large amount of data on this 
type of stirrer can be found in the literature. 
Wuirte and his collaborators found the following 
experimental correlation [18]. 


a ~ =" ¢ " (=) | ny 


MILLER and Mann [12] found the following 
relations : 


pasa’ 
‘_—_ 
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Fic. 16. Paddle stirrer. Influence of the ratio stirrer te 


vessel diameter on the dimensionless power consumption. 
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nas (“ 
fa v 


D cases. 
From the survey given by Hooker [7] it 
(in which only n and H were varied). This  fojjows that 
expression holds true for paddle stirrers with a 
blade angle of « = 90° as well as for paddle P (=) (¢ (5 


stirrers with a blade angle of «= 45°. The onde D D 


P i (7) proportionality constant is different in both 


v / 
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Fic. 17. Propeller stirrer and paddle stirrer with inclined 
blades. Influence of Reynolds number on dimensionless 


power consumption. 
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From our own investigations we found the ex- 
pression 


P nd? 0-15 d -1-38 
fe" (0) 


The influence of * and | was not further ascer- 


tained as these quantities have already been 
sufficiently studied. Hence, the following formula 
may be used with a fair amount of reliability for 
the power consumption of paddle stirrers : 


P nd2\-°™5 (d\-3 (H\* /w\?* 
‘ ~ 7 - 3 
pntd® ( y ) (5) (5 (7) ” 

These formulae only apply in the turbulent region, 

i.e. for Re > 10° to 10°. 


(Figs. 16, 17) 


For Re < 10? to 10° the flow becomes laminar, 
as was already found in the measurements of the 
mixing time. In this region 


P (= = 
pn®d5 v 
Turbomizxers—The relations holding for paddle 
stirrers also apply in the case of turbomixers. 
For turbomixers Mack [11] found the follow 


ing complicated relation in the turbulent region 
g Pp £ 


, 0108 


f 1.33 
F = orees(2) 
pn®d> d 


For practical use the formula may be simplified 
to 
\1 
Ba = 9-7 3° (=) (s = number of vanes) 
pn*d d 
The influence of the Reynolds number was negligible 
owing to the presence of stator rings. 

To ascertain what happens to the power con- 
sumption when the paddle stirrer is changed into 
a turbomixer the influence of vane angle 8 on 
power consumption was measured. The relation 
can be represented by the formula 

i Q)028 
pnéd® ~ (sin 8) 
This relation was ascertained for the turbulent 
region. 


Propeller stirrers and paddle stirrers with inclined 
blades (x < 90°). The formula (3) is also more or 
less applicable to propeller stirrers. 

Measurements by Sroops and Love. [16] 
produced the following result : 


P_ (nd*\°*/;d\°* (in these measurements 
% (;, only n and d were varied) 


pntd® \'y 
The diameter remaining the same, the power 
consumption increases about proportionately to 
the pitch of the propeller [13]. Only few reliable 
measurements on this point have been recorded. 
For paddle stirrers the influence of the blade angle 
a on the power consumption was ascertained (Fig. 
18). It was found that 
= ~ (sin gn 


pn*d> 


x 3- BLADE LE STIRRER 
@25x0-*m 


O 4- BLADE PROPELLER ACCORDING TO 
HIXSON AND BAUM 5) 


new 04 
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Fic. 18. Paddle stirrer with inclined blades. Influence of 
blade angle « on the dimensionless power consumption. 
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The relation between power consumption and 
Reynolds number is shown in Fig. 17. It is in 
good agreement with that for propellers, as given 
by Rusnton [18]. No baffles were used in any 
of these experiments. 

Baffling, eccentricity of stirrer—If the mixing 
vessel is fitted with baffles (see Fig. 14). or if a 
turbomixer has a stator ring round the turbine 


x duwe 6 xt2n10 4m? 
@ dxws 6x 1.2u 00! m2 Cece) 
DO duws 152 206-0 “‘nt ( 4 BAFFLES) 


+ 12-BLADE TURBO - mixer (f = 20, > 





wheel the influence of Reynolds number on the 
power number is almost nil. The power con- 
sumption increases by a factor of from 2 to 5 
(Fig. 19) when using baffles. 

Eccentricity of the stirrer has about the same 
effect on the power consumption as has baffling. 
Measurements on paddle stirrers showed an 
increase in power consumption by a factor of 



























































































































































Fic. 19. 








Paddle stirrer and turbomixer. Influence of 


Reynolds number on dimensionless power consumption. 
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from 1-5 to 2-5, this increase being 1 to 1-5 in 
the case of propeller stirrers (e = 0-5). 


3.6. Mixing energy 


Multiplication of the mixing time by the power 

consumption of the stirrer vields the energy 

consumed for complete mixing. 

a c 
P, Vy 
PPP Hy 
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SITUATION 1 
Fic. D 


SITUATION 2 


The theoretical minimum of mixing energy 
required can be calculated as follows (see figure). 
In situation 1 the potential energy of the two 
liquid layers with respect to the bottom of the 
vessel is : 


Evo, = Vi pif (he + By) + Ve pod (hg). 


In situation 2 the potential energy is equal to : 


poth 


+ 


; V, + pel, 
Exyotg = (Vy + V)( 1 P2" 2) g. (hy + hy). 
potg = (hy 2)( V,+V- g. Uh, 2) 
The increase in potential energy becomes : 


Ect =4 V; he (p2 — pie 


= ¢ E = 3¢(1 — ¢) Vh(p, — p,)8- 


The increase E represents the theoretical 
minimum of mixing energy required. It can reach 
a maximum of 0-125 Vh (pg — p,) g. In that case 
¢ = }, i.e. the two liquid volumes V, and V, are 
equal to each other. In our experiments this is 
the normal case. We found that this theoretical 
mixing energy is only a small percentage of the 
mixing energy actually required. 


4. DiscussIOoN OF 
OptTimuM STIRRING CONDITIONS 
4.1. Conditions for optimum performance (efficiency) 


With the correlations found for mixing time and 
power consumption optimum conditions can be 


derived for the dimensions and shape of stirrer 
and vessel. 
The optimum stirring conditions may, among 

other things, be defined as : 

(a) those in which the power consumption is 
lowest for a specified mixing time ; 

(6) those in which the mixing time is shortest for 
a specified power consumption. 


The data will be analysed according to the first 
definition. However, we would add that an ana- 
lysis based on the second definition in most cases 
leads to the same results. 

For a given mixing process the physical pro- 
perties of the liquid are, of course, fixed. The 
shape and dimensions of stirrer and vessel are 
generally, within certain limits. left to the 
When these variables, too, have been 
established, the speed of rotation must be such 
as will ensure mixing within the required time. 

For the mixing time T we may write : 


Ee =; 
» Vz , , 9 
pa” 8) 


designer. 


T = T (n, 


and for the power consumption P : 


P P (n, dye 


D pa?) 


We will successively deal with the influence of 
the volume of the mixing vessel V ; 


; ¢ . a 
the stirrer diameter ratio D' 


, ; —r H 
the height ratio of the liquid in the vessel D' 


the width ratio of the stirrer blade ft 


blade angle a; 
vane angle 8. 


The procedure for determining the optimum 
conditions consists in varying n and only one of 
the other variables, the remaining ones being kept 
constant. For a constant mixing time the variation 
of n depends on that of the variable that has not 
been kept constant. In this way the value at 
which power consumption is lowest can be deter- 
mined for any of the variables involved : Because 
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of the special structure of the formulae for T 
and P the optimum value found for one variable 
is little influenced by even quite considerable 
changes of the other variables. 

4.2. Influence of liquid volume V in mixing vesse 
pn*d? 0-3 


A pg 


If + ~| and, moreover, 7 does not 
depend on the Reynolds’ number equation (2 


may be T ~n"'(n*d)°*, while 
0.15 


written as 
equation (3) becomes P ~ n*d5 (nd*) 


Elimination of n gives : 
P~T 8 d*"" or 
P — T 1.78 i 4 


(as d* ~ V) 


Power consumption per unit volume, P/V, is 
proportional to T~'™* V°* and therefore increases 


when the liquid volume increases (no optimum). 
, 


The mixing energy per unit volume, » which 


is proportional to T° V°* also increases at 
higher V values as well as at smaller 7 values. 

pn*d? 
Ap gH 
r, P~T?* V'* and 


ye, 


When the quantity has no influence on 


> 


therefore w Te 


It is therefore advantageous to use a small 


mixing vessel. In large vessels several small 


stirrers pumping in the same direction are to be 
preferred to one large stirrer. 


4.3. Influence of stirrer and vessel dimensions 


Influence of the stirrer /vessel diameter ratio 


; (*, and DPH are constant ) 
The equation for T becomes : 


d*® P 
pa) (Re > 10°) 


T~ nregres| 
a 
where f (ip 5 783) is given by Fig. 9, and 


equation (3) becomes 


P~ nes aes (fe 


f da’ ze . 

For pe 7 < 0-04 I (pes = = constant. 
For the mixing time T it is then true that 
T~n'**a**, Furthermore, P ~ n®* d*'. 
Elimination of n produces P ~ T*'* d**. For 
a given time 7 the power consumption will 
decrease as the stirring diameter increases. Hence, 
the best plan will be to take the largest possible 

3 
For prs yet 
d® d* 
I ( pes qs) ~ Fe 


becomes : 


T~nted**| 


diameter. >» 0-08 it is true that 


The mixing time T 


d 8 
4 or T~n'*da?* 


. i 
] 2.85 i** 
P~n r (|, 


and 
Elimination of n produces : 
13 
P~ r+ gs (S) o P~T'" &*. 
D 


In this region P will increase with increasing 
stirrer diameter if T is kept constant. The best 
plan is then to keep the stirrer diameter as small 
as possible. 

There is, therefore, an optimum diameter in 
the region : 


d3 d\ H 1/6 
0-04 < p37 * 0-08 or 0-34 < (5)( p) < 0-43. 


Practically speaking, this means that this 


, d 
diameter corresponds with a value of — & 0-4, as 


D 1/6 
the value of (|) in most cases will not differ 


very much from 1. 

The above results are only valid for a turbulent 
field of flow. In the laminar range (Re < 250) the 
mixing time very rapidly increases at lower 
Reynolds’ values (Figs. 7 and 8). In this case 
it is better to take the largest possible stirrer 
diameter, as will be shown below. In this region 
it has been assumed that : 

Tnd* { pn*d* \°* = /pnd* f | €& 

(Sy ~() shee} 


and : 
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we~("S) -() 


If a ) “iY f ¢ > 0-08 and 
Sea (p) ( ” pm) > 0-08 an 
H 


= constant], we find : 


P~T°* ¢"", 


which proves that here a large diameter is 
favourable. In actual practice the stirrer dia- 
meter will have to be about equal to the vessel 
diameter when mixing viscous liquids, 

4.3.2. 


Influence of liquid height = (5 and D*H 


are constant } 
In the same way as shown in the preceding section 


H . ; 
the influence of pn be found. From equation 


2y7)j2\ 03 / 
(2) it follows that T ~n'd* (“ a £(5) 


while from equation (3) it is found that 


0-6 
P ~ n*d° (nd*)°™ (5) . Elimination of n now 


(DH is constant) gives the following relation : 


ve przez! ¢( ‘o d® - 0 
P~T*H (n)} For pps qos < 0-05, 


s(7) = constant, so that then P ~ T' H', 


In other words, P increases with increasing H 
values. 


7 d* . H i\-** 
For D?* 7s > 0-05, t(5) (7) and as 
D~ HH, P~T 18 H'? 1-35 T 1-8 Ho}, 
Hence, in this region P slightly decreases with 


increasing H values. There is therefore an 
3 


, ; a ; nn Os 
optimum value for p®* 7 = 0-05. When D is 


H 
taken equal to 0-4 the value for p Preves to 


become * 0-8, which is a normal charge depth in 
actual practice. 


4.3.3. Width of stirrer blade (influence of “) 


From the formula for power consumption and 


mixing time, giving the influence of the width 
of the stirrer, we may derive rules for the most 
favourable dimensions of the stirrer. 


F ° . P w\** 
rom eq. (3) it can be derived that ond ~ (5) , 


while T ~ (nd*w)" holds good (if = < 0-5), when 


the Froude number has no influence. Elimination 
of n gives P~ T*w**®. Hence, P ~ w*® for 
T = constant. 

According to this formula the blade should be 
made as wide as possible. 


w w 0-66 
For values of F > 0-5, T ~(nd*wy" (") 
( 


(see section 3.2.1), hence, T ~ n™' w For 
the power consumption P we find after elimina- 
tion of n: P~ T* w®, 


Thus, for 


0-33 


5 > 0-5 there is still a small effect 
a 


w - , 
of — on the efficiency. In practice the ratio 


( 
w 
d 
When the density Froude number exerts an 


= 0-5 is never exceeded. 


influence we obtain for 


© 05 PY TM a 


>05 Pw~T'* w 


d 


i.e., the trend remains the same. 


4.4. 
Influence of blade angle « 


Type of stirrer 
4.4.1. 


The value of the blade angle « more or less 
determines the type of stirrer (paddle or propeller) 
and the influence of « shows which type of stirrer 
is the most effective. From the relation found 
between mixing time and blade angle it can be 
ascertained what blade angle is most suitable. 

T ~n"*(sina)™ k, =1to2 

P ~n?* (sina) ky = 2to8 
Elimination of n produces 


P ~ T* (sin a) 97" 
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For constant mixing time T, P varies according 
to P ~ (sin x)? "*"_ For (ky — 1-7k,) > 0, P 
with for (kg — 17k.) < 0 it 
From Figs. 14 and 18, ky 
and k, for various values of « may be found. 

For 0 <2< 45, kh, ® 
that A, — 17h, ~ — 02. 
k, = 2 and k, = 2-4; then ky, — 17h, ~ 0. 

The small values found for k, — 1:7 k, make it 
clear that the influence of « is only very small; in 
other words, the propeller stirrer has about the 


increases sin a; 


decreases with sin =. 


1 and ky > 1:5. so 
For 45° < « < 90°, 


same performance as the paddle stirrer. 


4.4.2. 
The influence of the vane angle decides whether 
a paddle (8 90°) is 
to be preferred. 


Influence of rane angle B 


90°) or a turbomixeir (8 - 
The equation for the mixing 


time in this case becomes : 
T ~n™'* (sin py" 
and the equation for the power consumption : 
P ~n'*™ (sin By. 
Therefore, P ~ T~'* (sin By. 
Hence, 


mixing; in other words, the paddle mixer is 


a high sin 8 value is favourable for 


more eflicient than the turbomixer. 


4.5. Comparison with jet mixing 


Jet mixing with external circulation of liquid is 
An important 
Fosserr and 


frequently applied in industry. 
point is the size of the nozzle. 
Prosser [4] investigated into the possibilities of 
jet mixing and derived an expression for the 
mixing time. This enables us to compare jet 
mixing with mixing by agitation by extrapolation 
of the of Fosserr and 
ProssER to our small-scale experiments. The 
basis of comparison is again power consumption 


large-scale results 


with equal mixing times (see Fig. 20). 





H NOZZLE 
VELOCITY Vv 











VAN DE VuSsSsE 


The velocity of the liquid v at the nozzle must 

v k H (py, - Pa) (4), 

2g py sin® @ 

Pa “ Pr 
P2 

0-079, k= 14, 


comply with the condition 


where & is a constant still dependent on 


f ~ (° *) , |. For Pa pi 
P2 P2 


H 118 mm. Let sin @ be 
1 


then 


by2, 
The mixing 
sp? 
V & 
0-120 Q°* sec (Q in m® /sec). 
The required input is found from 
P=tpuv®x@Q 10°*QNm sec. Fig. 20 
gives the P and T values corresponding to various 
@ values. It will be seen that the power con- 
sumption in jet mixing is about four times as high 
as in paddle mixing. With the aid of a P, V vs. 
T graph (Fig. 20) an adequate comparison between 


026m and v 2-30 m. sec. 


vo 


can be found from T For 


time T 
dD 148mm. 7 


power 


2-6 
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Fic. 20. Comparison between jet mixing and paddle 


stirrer. 
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various stirring devices can be reached. Elimina- 
tion of Q gives PT? = 37-5 Nm sec. 


4.6. Concluding remarks 
A general survey of the efficiency of the various 
types of stirrers, using the optimum. stirrer 
dimensions, is given in the table below. 


5. Continuous MixinG or LiquIpDs 
5.1. Introduction 


Continuous mixing of liquids is understood to be 
the operation where one or more streams of liquid 
are continuously fed into and discharged from a 
mixing vessel. Both the flow rates and the com- 
positions of the liquids may vary with time. The 
flow rate at the outlet will at any time be equal to 
the sum of the flow rates of the inflowing streams ; 
the final composition will be a function of the 
flow rates and compositions of the inflowing 
liquids and of the degree of mixing in the mixing 
vessel. 

Mixing in the mixing vessel is called perfect 
(one ideal mixing stage) if the inflowing liquids are 
immediately and completely mixed with the fluid 
in the vessel, so that its properties are uniform 
and identical with those of the outgoing stream. 
This case is frequently encountered in the 
literature and calculations on stirred reactors are 
often based on this assumption. 

Complete absence of mixing occurs if the 
elements of the entering liquid move through the 
vessel at constant and equal velocities along 
parallel paths, so that no mixing takes place 
either in the direction of flow or in other (radial) 
directions. 


The case of no mixing in the directions of flow 
(perfect piston flow) can be approached by an 
infinite number of ideal mixing stages lined up 
in series. 

The normal case of a non-perfect mixer is very 
difficult to handle mathematically. The concep- 
tion of an eddy diffusivity in axial and (or) in 
radial direction as usually applied in packed beds 
and in particulate fluidized beds, does not hold 
for the case of an agitated vessel. 

The continuous-flow mixing of one liquid 
stream can be suitably characterized by the curve 
of the outlet concentration against time, if the 
inlet concentration is a well-defined function of 
time (e.g. a step function). 


In many processes two liquid streams with 
constant concentrations must be mixed. In this 
particular case only mixing in a radial direction 
is important, which can be effected, for instance, 
by line blending. As a rule, the liquid will be 
uniform after having passed through a pipe length 
of about thirty times the pipe diameter (turbulent 
flow conditions). Bends and constrictions appre- 
ciably improve mixing. 

To conclude our investigation, a few experiments 
were performed with one liquid stream to evaluate 
the relation between batch mixing and continuous 
mixing. 


5.2. Measurements (mixing times) 
The experiments were carried out in a mixing 
vessel (H = D = 0-15 m) with a centrally placed 
stirrer (Fig. 214). In order to obtain an even 
distribution of the inflowing liquid stream the 
vessel was provided with an inlet compartment 





Comparison of some types of stirrer on a basis of power consumption (in arbitrary 
units) required to ensure a certain mixing time T Optimum dimensions of stirrer 





Paddle stirrer 
Turbomixer 
Propeller stirrer 


Jet mixing 
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consisting of perforated plates with the same 
diameter as the vessel. In this way mixing caused 
by the injection of the fluids was suppressed. 
Disturbance of the flow pattern by the outlet 
was avoided by using a sieve bottom. 

As the sole object of this investigation was the 
study of stirrers, the influence of different types 
of inlet nozzles, ete., was not further examined. 

A small amount of dye solution (bromophenol 
blue) into the 
inflowing liquid stream (10 ml in 1 sec). The 
distribution of the colour through the mixing 


was instantaneously injected 


vessel and in the outgoing streams (under the 
The outlet 
and a 


sieve bottoms) was watched. con- 


centration curves for an ideal mixer 
completely non-ideal mixer are given in Fig. 218 
(curves 1 and 2). 

The outlet concentration will show a maximum, 
which in the normal case of a non-perfect mixer 
will appear a time 7 after the moment of dye 
injection, which time is shorter than the residence 
ti (oie of = (Fig. 21m, curve 3). 

flow rate 

In most experiments an even colour distribution 
throughout the vessel could be observed some 


time after the injection (2), which “* mixing time ” 
Fic au 


PERFORATED INLET 
COMPARTMENT 
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Fic. 21. Continuous mixing. 


corresponds to the time T. The piston flow mixer 
(Fig. 218, curve 2) will, of course, never give such 
a “ mixing time.” 

The mixing measured as a 
function of agitator speed and flow rate. The 
influence of the flow rate on T in the case of a 


times 7’ were 


stagnant stirrer is shown in Fig. 22. The times 


then found (called 7) are always shorter than 


i= 
, will, 
R 


of course, depend on the conditions of flow in the 


the residence times (7',). The value of 


mixing vessel, hence, on the construction of the 
inlet, channeling effects, ete. 

The influence of the stirrer speed on the mixing 
time 7’ when there is no flow is the same as in 
batch mixing (mixing time 7’,). 

In the case of both liquid flow and stirring the 
mixing times 7° appear to correlate with the 
batch mixing time 7, and the time T,, as is 
The best correlation for T i 
Re ae 
. 


shown in Fig. 23. 

terms of T, and T, is given by 
0-5 

VT; T: 
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Fic. 22. Continuous mixing. (Mixing) time Ty, vs. flow 
rate Q (no stirring). 
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Fic. 23. Continuous mixing. 


Influence of the batch mixing time 7, (only due to stirring) 


and the (mixing) time T, (only due to the flow) on the mixing time T (stirring and flow). 


Summarizing our results, it may be said that in 
the case of continuous flow through a stirred vessel 
the mixing time T can be predicted from the mixing 
time in batch operation (7) and a time T,, 
which is always shorter than the residence time. 

If the batch mixing time 7, < Tp, then 
T <Tp, in which case a perfect mixer is 
approached. 
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es = 


NOTATION 


stirrer diameter 


= vessel diameter 


eccentricity of stirrer 

acceleration of gravity 

height of liquid in vessel 

number of revolutions of stirrer per unit 
time 


= refractive index 
= pitch of propeller stirrer 


power consumption 
pumping capacity 

radius 

number of vanes (or blades) 
mixing time 


= velocity 


liquid volume in mixing vessel 
width of stirrer blade 
blade angle (paddle with inclined blades) 
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B = vane angle (turbomixer) 
p = dynamic viscosity of liquid 
vy = kinematic viscosity of liquid 
p = density of liquid 
w = angular velocity 
TQth 
— 
pn? d? 
Apgll 
pnd* 
mv 
a 
pPs °° 


dimensionless mixing time 
density Froude number 
Reynolds number 


= dimensionless group comparing the 
dimensions of the stirrer and the dimen- 
sions of the mixing vessel 


pn? dd dimensionless power consumption 
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Uni-directional heat and mass transfer at high humidities in a 
wetted wall column 


R. C. Carrns and G. H. Roprr* 


School of Chemical Engineering, N.S.W. University of Technology, High Street, Kensington, New South Wales, 
Australia 


( Received 9 January 1955) 


Summary—Further work is presented on simultaneous and adiabatic heat and mass transfer 
at high humidities. In previous work [1] the authors studied heat and mass transfer in opposite 
directions in a wetted wall tower. In the present study the authors have investigated the 
transfer of heat and mass in the same direction, in a similar wetted wall column. 


Data were obtained for the range 0-49 to 0-80 for the mol fraction of the diffusing vapour 
in the gas film. 


The equations of Cutron and Cotsurn [2] and Co_.surn and Drew [4] which were found 
to be not applicable for adiabatic humidification at high humidities, are also not applicable to 
the new data. The mass transfer data could not be correlated as a function of (pyyy, P) as 
before, since the effect of a, the ratio of the nett sensible heat carried by the diffusing vapour 


to the heat transferred in the absence of mass transfer, had to be included. 
The heat transfer data for both adiabatic humidification and dehumidification were correlated 
by 


1 on 
(h/hy) {1 + >) = (ppy/P)** 
which may be written 


(h/eG) > Pr? ® (pyyy/ PY?) (1 4 < 0-025 Re? 


Résumé — Les auteurs présentent un nouveau travail sur le transfert simultané et adiabatique 
de chaleur et de masse pour des humidités élevées. L’étude se fait dans une colonne a parois 
mouillées comme précédemment, mais avec les transferts de masse et de chaleur s’effectuant dans 


le méme sens au lieu de s‘effectuer dans des directions opposées. 


La fraction moléculaire de la vapeur diffusant dans le film gazeux correspond a Vintervalle 
0,49-0,80. 


Les équations de Curiron et Cotpurn et de Cotpurn et Drew qui ne s’appliquaient pas 
4 humidification adiabatique pour des humidités élevées, ne s‘appliquent pas d'avantage aux 
nouvelles déterminations. On ne peut pas, comme précédemment, représenter les valeurs du 
transfert de masse par une fonction de (pyyz/P) puis-qu'il faut tenir compte du facteur a,” 
rapport de la chaleur sensible, nette apportée par la vapeur diffusante 4 la chaleur transmise en 
absence du transfert de masse. 


Pour une humidification et une déshumiditfication adiabatiques, le transfert de chaleur peut 
étre représenté par les deux équations ci-dessous : 


(h/hy) (1 + 3) =(pam/P?™ 


(h/cG) - Pr? (ppy,/ PY? ™, (1 + 3) = 0-025 Re®? 
* Present address : Department of Chemical Engineering, University of Delaware, Newark, Delaware, U.S.A. 
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INTRODUCTION 


Tuts paper is a continuation of previous work [1], 
in which the authors investigated the adiabatic 
humidification of air at high humidities in a 
wetted wall column. Heat and mass were trans- 
ferred simultaneously in opposite directions with 
countercurrent flow of the liquid and gas streams. 
Previous correlations of GILLILAND and SHER- 
woop [5] and Cuitron and CoLBuRN [2] did not 
apply at high humidities. 


For mol fractions of water vapour from 0-08 to 
0-85 in the gas film, correlations by new heat and 
mass transfer factors gave : 


Jn = (h/eG)- Pr’? (pang : aaa 
ip = (ke Paw M,,/G) - Sc? - (ppy/PY?" 
0-025 Re®? (1) 
for a liquor rate of 77 lbs/hour in a 0-90 in. in- 
ternal diameter column. 


Co_purN and Drew [4] theoretically developed 
the relation 


q.=h At 


a 
aa hy A At (2) 


to predict the effect of mass transfer on the heat 
transfer coefficient. A, is the heat transfer co- 
the absence of simultaneous mass 
is the 


efficient in 
transfer. The dimensionless number, a, 
ratio between the nett sensible heat carried by 
the diffusing vapours and the heat transferred in 
the absence of mass transfer. This is negative if 
mass and heat are transferred in opposite direc- 
tions and positive if mass and heat are transferred 
in the same directions. 

The authors have previously shown that their 
data for heat transfer could not be correlated in 
the form of equation (2). 

In order to test the correlations previously 
developed by extending the value of a to the 
positive region, it was decided to operate the 
wetted wall column so that the gas stream was 
simultaneously cooled and dehumidified. The 
theory of CoLpuRN and Drew predicts that the 
heat transfer coefficient will be increased by 
vapour diffusing in the same direction as the heat 
transfer. 


EXPERIMENTAL PROCEDURE 


AND RESULTS 


The original apparatus was modified slightly to 
permit the study of unidirectional heat and mass 
transfer. The make-up water system was replaced 
by a vessel provided with an overflow to main- 
tain a constant amount of water in the liquid 
recirculation system. The water vapour con- 
densed in the column was measured by the 
amount drawn off from this overflow. Apart 
from the measurement of the overflow in place 
of make up, the procedure was identical with that 
previously described except that the entering 
water temperature at the top was maintained 
lower than that leaving the bottom. An inherent 
feature of the operation of this apparatus is that 
for mass to be transferred to the water simul- 
taneously with heat, the bottom and top tem- 
peratures cannot be matched. 

Further modifications included the installation 
of a larger condenser and the 0-9 in. column 
previously used was replaced by a tube 1-04 in. 
internal diameter and 37} in. long. A gear pump 
was used on the liquid recirculation system in 
conjunction with a constant head feed in place 
of the centrifugal pump, which initially gave 
trouble due to leaks in this work. 

A calibrated Smith dry gas meter, 0-100 c.f.h., 
was used to measure the volume of atmospheric 
air entering the column together with a rotameter 
for control. The water system flowrator was 
calibrated for temperature variations since the 
sensible heat gained by the water is, in this case, 
equal to the enthalpy change of the air stream, 
and the water rate is required to be known 
accurately. 

The experimental results obtained are shown 
in Table 1. 


METHODS OF CALCULATION 


The runs were interpreted in the same manner as 
in the previous work. The mass transfer rate was 
taken as the mean between the overflow and the 
difference between the total water vapour entering 
and leaving the column. The total water vapour 
entering the column is the sum of the water 
vapour from the boiler and the incoming air 
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Table 1. 
| 


Experimental results. 





Il air | Boiler 
lbs /lb dry air 
rate - - 
(lbs /- 


hour) 


Dry 


air 
| 


Run 
No. 


rate 


In Out 


4-238 OWT “O135 5-70 "148 


3-055 ‘OO O168 "180 | 4°127 


O1LST “234 4-024 


“O106 175 4-975 


OSLO “150 4564 


“0260 ‘191 | 5-104 


3-982 03543 O95 | 5297 


6-420 O40 4-689 


0223 


6-354 | “038 5-018 


6-373 “086 4491 


| 6-334 | 0305 5-083 


3-894 0225 4-649 


1-471 “0168 4-708 


1-466 | 0005 “0153 3-798 | 


0097 | 0190 | 


| 3-684 4-000 | 


3311 | 0006 | -o182 | 3-995 


Con- re Over- | Mixture Temp. | 
rate densate | rate | flow 
(lbs /hr) | (lbs /hr) \(Ubs /hr) | (bs hr) 
‘ —+—_- 
4°283 | 1-213 


1-242 


‘716 


‘734 


1-605 
1-463 | 


1-500 


rate —_— — 
| 


Column) Liquor | 


rate at 
top | 


Water temp. 
at column 


at column °F | pres. 





Top | Bottom atm. Bottom 
- | _ 





—_|___ 
266 348 1015 | 73: . 187-4 





318 | 1-013 | 


73-9 192-0 


196-7 


262 
261 10138 | 74-0 
| 





| 1000 | 740 186-9 196-7 
| 





| 


1007 | 71-0 175-3 190-2 


203-0 


193-9 


1-004 | 


1-004 71-2 180-6 


70-9 





156-3 181-9 


164°8 


1-007 


1-008 710 | 184°3 


221 261 1-007 708 160-8 178-7 


71-0 


256 1-007 176-6 184°5 


221 


265 


257 


226 1-005 190-9 





230 1-004 200-6 





228 257 | 1-004 199-1 
| 


230 1-004 


226 “202. | 1-008 | 71-1 | 166-1 
| 





272 186-4 





189-0 








- All calculations were carried through to four significant figures to avoid cumulative errors associated with the 


length of the calculations. 


the steam 


line. The total water vapour leaving the column 


corrected for the condensate from 
is the sum of the water vapour in the exit air and 
the water from the condensed product. The mass 
balances are expressed as an error in closure 
defined by : 


total vapour in —total vapour out 


%, error in closure = x 100 


overflow 

The inlet humidity was calculated from the dry 
air rate and the mean of the total water vapour 
entering the column and the total water vapour 


leaving the column plus the mean mass transfer 
rate. The outlet humidity was calculated from 
the dry air rate and the mean of the total water 
vapour leaving the column and the difference 
between the total water vapour entering the 
column and the mean mass transfer rate. 


The driving force for mass transfer was taken 
as the logarithmic mean between the partial 
pressure driving forces at the top and bottom of 
the column. In order to evaluate the true mean 
driving force it would be necessary to find the 
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operating line for the cooling process. COLBURN 
[3] suggests using the equation : 


dp, _(P — p,) (Pp — Pe) (€* — 1)- Pr” 
dt, 7 Ps (t, = t.)a §,.24 


This equation was derived on the basis of the 
CutLtton-CoLBuRn [2] analogy and the CoLBuRN- 
Drew equation (2). In adiabatic humidification 
at high humidites it has been shown that these 
relationships do not explain the data. Hence 
there is little reason to expect that these relation- 
ships would hold for dehumidification. As will 
be seen later, the results of this investigation 
cannot be correlated by the CHtILToN-CoLBURN 


(3) 


Table 2. Calculated 


equation or the CoLpurnx-Drew equation and 
that a general relationship between heat and 
mass transfer coefficients at high humidities is 
not known. Hence true mean driving forces 
cannot be calculated. The use of the logarithmic 
mean driving force was based on the close equality 
of this mean and the true mean for both partial 
pressure and temperature during adiabatic 
humidification [1]. 

The calculated mass transfer data are given in 
Table 2. 

The heat balances were calculated from a 
datum of 100°F. The water rate at the bottom 
of the column wes equal to the sum of the rate 


Mass Transfer Data. 





| H at column 
Run lb /b dry air 


Ne. |-——_,—_ —_- 


Ny | 
(lb moles 
_howr) 


(AP) mn 
(atm.) Re 



































11 328 | 





a | vane) 





13 | e210 | 


14 | 2-615 | 








1s | 1-110 | 





(kg RTd) 
an 


O1l4 


Pam (*cPem)\ .2s 2 
Re oes Po (grat) se RO 


O177 






























































| 
16 1-227 | 
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at the top of the column and the mean mass 
transfer rate. The heat balances are expressed 
as an error in closure defined by : 


*% error in closure = 
— heat change in the water — 
enthalpy change in the air stream x 100 
enthalpy change in air stream 
In calculating the convective heat transfer 
coefficient the heat transferred by convection 
was evaluated as the sum of the sensible heat lost 
by the dry air, the sensible heat lost by the total 
water vapour entering, in cooling from the inlet 
to the outlet temperatures, and the sensible heat 
lost by the mean mass transferred in cooling from 


the exit gas temperature to the average liquid 
temperature. This quantity was corrected for 
radiation as before and logarithmic mean tem- 
perature differences were used. 

The calculated heat transfer data are shown in 
Table 3. 


Discussion 
The results were plotted by several of the methods 
suggested by previous investigators. The liquor 
rate varied only slightly and the correction to the 
standard liquor rate of 77 lbs /hour was negligibly 
small [1] (Fig. 2). GiLLiLanp and SHEeRwoop [5] 
correlated their results at low humidities by a 
Nusselt type equation which for a Schmidt 


Table 3. Calculated heat transfer data. 





qT (Add im. 
(Btu/hr) | ( Btu/hr) (°F) 


h 


350 


328 126 














) PP® RO? 


Heat 
balance %, 


0231 


0210 








0199 
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number of 0-55 and a liquor rate of 77 lbs/hour 
may be reduced to 


(kg RTd/D)(ppy/P)- Re’ = 0-017 (4) 


The present authors correlated their data for 
adiabatic humidification, for high humidities at a 
liquor rate of 77 lbs hour by the Nusselt type 


equation 


(kg RTd/D)(ppy/P)™ > Re’™ = 0-016 (5) 


The results for dehumidification are plotted in 
Fig. 1 
below the lines of equations (4) and (5). 


seen to fall consistently 
The data 
are too scattered to permit a correlation in the 


and the data are 


form of equations (4) and (5) 


0-040 





‘.| GILLILANO 
+. , SHERWOOD 
‘¥ 


ADIABATIC 

HUMIDIFICATION 

= ne 
ie 


| 
, 


“St 


Fic. 1. Comparison of the mass transfer data with equations 
(4) and (5). 











Cuitton and CoLpurn [2] correlated the heat 
and mass transfer data at low humidities inside 
tubes by equations which at a liquor rate of 
77 Ibs “hour are 


(h cG) Pr . (ke Papo M,, G) ‘ Sc* , 


0-026 Re®* (6) 


Cairns and Roper [1] correlated their heat and 
mass transfer data by equation (1) for high 
humidities. The experimental mass transfer data 
for unidirectional transfer, as shown in Fig. 2, 
fall below the lines of equation (1) and (6). The 
points in Fig. 2 are scattered and a satisfactory 
correlation of the data in the form of equations 


. Carans and G. 


H. Rorer 





0-030 


scar COLBURN 


. ADIABATIC 
HUMIDIFICATION 











010 , Pen 
, 


100 


Fic. 2. Comparison of the mass transfer data with equations 
(1) and (6). 


(1) and (6) is not possible. It was assumed that 


there is a relationship between 
- Re®? 
and log (pyy/P) + ba 


log (ke Peo M,, - Sc? . 
The data were best fitted 
The plot of these 
groups is shown in Fig. 3. It was not possible 
the data of adiabatic 


where 6 is a constant. 
when 6 was taken as 1°85. 
to correlate humidification 
and dehumidification by this type of correlation. 
The 
line shown in Fig. 3. 
in Fig. 3 is the point representing the CurLron- 


humidification data lie above the 
The theoretical point shown 


adiabatic 


CoLBURN correlation for the case of dilute mix- 
tures, i.e. when (ppy/P) = 1, anda = 0. 

The heat transfer data have been plotted in 
Fig. 4 in the form of (h/cG). Pr*/® Re®* versus 
(Ppm/P) as suggested by Carrns and Rover. The 
data fall consistently below the Cu1LToN-CoLBURN 
line and above the line representing the data for 
adiabatic humidification. No correlation was made 
of the data in this form but it is evident that 
higher coefficients are obtainable with positive 
values of a than for negative values, as in adia- 
batic humidification. 
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; 





















































-O%8 -0O4 
LOG ( "t+ 1880 
Fic. 8. Mass transfer data plotted as a function of 


log (Pam + 1-85a). 


-02 


The ratio h/hy versus a is plotted in Fig. 5. As 
was the case for adiabatic humidification, it was 
not possible to correlate the data. The theoretical 
line of CoLBURN and Drew is shown also in Fig. 5. 
The Co_purn and Drew equation can be re- 


arranged to give 
a a 


h/ho = ——~ + 1+ (7) 


5 2 
for the values of a used in this investigation. In 


Fig. 6 the dimensionless group (h /hg) ! ‘) has 





0 0b ~ 











a 
CHILTON - cCoLeur 


b  ADIABATIC 
HUMIDIFICAT 

















Fic. 4. Comparison of the heat transfer data with equations 


(1) and (6). 





T T 
COLBURN - DREW 
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%> ioe) 02 O4 
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Fic. 5. Comparison of the heat transfer data with the 
theoretical equation (2) of Co-puRN and Drew. 


been graphed against (pyy/P) for both adiabatic 
humidification and dehumidification. The whole 
of the data can be well correlated by the equation 


h/he = (1 


which can be written 


+ ;) (Pam/P)™ (8) 


(h cG) ? Pr*’® - (pang Pye? (1 = 2 5) 
= 0-025 Re®? (9) 














ee 
COLBURN DR 
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THIS INVESTIGATION 














| 
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fae = ADIABATIC “HUMIDIFICATION _ 
- , DEHUMIDIFICATION a 


0-10 Pan 
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1-00 


Fic. 6. Comparison of the heat transfer data for both 
adiabatic humidification and dehumidification with 
(Ppm/P). 
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CONCLUSION 


Mass transfer coefficients for the case of adiabatic 
dehumidification can be correlated by inclusion of 
the groups (p,y,/P) and a and are not a function 
of (ppy/P) alone as was the case for adiabatic 


humidification. The presence of heat transfer 
appears to have reduced the mass transfer rate. 
The presence of mass transfer, when occurring 
in the same direction as heat transfer, increases 
the heat transfer rate, as was predicted by 
CoLBuRN and Drew. These increased heat trans- 
fer coefficients can be correlated only by inclusion 
of both groups (pyy P) and‘a, and are not a 
function of a alone as was derived by CoLBuRN 
and Drew. 
humidification and dehumidification. 


The correlation is made for both 


NOTATION 

the ratio between the nett sensible heat carried 
by the diffusing vapours and the heat trans- 
ferred in the absence of mass transfer. 
the area of the gas-liquid interface, taken as the 
inside area of the column, when dry, sq. ft. 
a constant. 

= specific heat, Btu /(lb.)(°F). 

= diameter of column, ft. 
molecular diffusivity of the diffusing component 
sq. ft. /hr. 

= base of Naperian logarithms. 

= mass velocity, Ib. /(hr.) (ft.*). 

= convective heat transfer coefficient, Btu/(hr.) 
(ft.?) (°F.). 


H. Roper 


ho 
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The solution of multicomponent distillation columns 
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Summary—The equations of A. J. V. Unperwoop for solution of distillation columns are 
arranged to give a solution in the form of a sum of difference products. For many distillation 
problems the present solution is much less tedious than the straightforward solution of 
UNDERWOOoD's equations or the determinantal solution proposed by ALpeR and Hanson [1]. 
The solution is illustrated by an example employing a five-component system. 


Résumé— Pour la solution des équations de A. J. V. UNpERwoop permettant le calcul des 
colonnes 4 distiller, les auteurs conseillent l'emploi d'une somme de produits d'une différence. 
Dans de nombreux problémes de distillation, cette solution est beaucoup moins pénible que 
la solution directe des équations d’ UNDERWOOD ou que la solution sous forme de déterminants 


proposée par ALDER et Hanson. 


La solution est illustrée par un exemple sur un systéme a cing composants. 


Tue solution of distillation 
problems is, in general, time-consuming and 


tedious work. The methods which are available 


multicomponent 


for such solutions employ various basic assump- 
tions and various techniques, and the particular 
method which is best suited to the problem at 
hand depends largely on the accuracy of the 
answer desired and on the physical make-up 
of the system to be distilled. As an example, 
plate-by-plate — calculation 
distinctly useful for some systems, for example 
systems in which all of the components of the 
feed appear in one or the other of the products in 
well-definable amounts and in which the difference 
in volatility between the key components is 
large. For such systems, calculation plate-by- 
plate from one end of the column will provide 
the solution within a few plates. If, on the other 
hand, the key components lie intermediate in 
volatility within a large number of components, 
and a large number of plates is to be expected, 
plate-by-plate calculations become difficult, and 
procedures which allow the calculation of the 
number of rectifying and stripping plates as 
complete sections become most useful. Such 
“group” procedures have only recently been 


procedures are 


developed, and to the authors’ knowledge, four 
basic solutions presently exist. The relative 
merits of these are as yet undetermined but one 
which shows promise of considerable usefulness 
is the set of equations derived by A. J. V. 
Unperwoop [8, 4, 5]. 

The equations provide exact solutions of 
distillation columns under the basic assumptions 
of constant relative volatilities and constant 
sectional flows. The correctness of these assump- 
tions depends upon the system, and although 
they are not in general exact, distillation equations 
involving their use have found wide application. 
The flows are usually set by the internal reflux 
at the top of the column and the thermal condition 
of the feed. Various procedures are in use for 
estimating average relative volatilities. Ex- 
perience has shown that satisfactory values may 
be obtained using an average of the top plate 
and reboiler temperatures weighted by the molar 
amounts of top and bottom products or, alter- 
natively, that temperature at which the product 
of the equilibrium constants of the light and 
heavy key is equal to unity. 

The solution of UNpERWoop’s equations gives 
directly the number of rectifying and stripping 
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plates. The systems ordinarily encountered in 
multicomponent distillation contain components 
both lighter than the light key and heavier than 
The 


following the lines proposed by 


the heavy key. trial and error solution 
, UNbDERWOoOD 
[5] then becomes prohibitively time consuming. 

Atper and Hanson [1] 
solution of the 


form, bypassing the calculation of feed plate 


have proposed a 
equations in determinantal 
concentrations proposed by UNpERwoop. How- 
ever, since the order of the determinants which 
the method of ALDER 
Hanson is equal to the number of components 


must be solved in and 
in the feed, the difliculty of the solution mounts 
rapidly with an increase in the number of com- 
ponents, and a solution which would avoid this 
difficulty is considerably more useful. The present 
paper proposes a solution which accomplishes 
this by developing the determinants into a 
sum of terms, each of which contains a quotient 
of products of differences, 

For the of the 
introduce the notation to be used, a brief review 


of the developments which led to the solution 


convenience reader, and to 


is presented. To avoid excessive duplication of 
previously published material, the review and 
the development of the solution will be presented 
for a specific case. The use of the solution in 
solving other cases can be obtained from a brief 
study of the previous paper by Alper and 
Hanson. 

For a system of three components A, B and C, 


then, the following equations may be written : 


apd (te), 


aco 


ai a 44 (@4)4 , %#f (7 p)a 


1 
aoe ame m0) 


a4b(r4)p 


by apd (2 p)p . ach (re), 
t4—¢ 


ay — aco 


— V’ (2) 


where V = rectifying section vapour flow, moles 
V’ = stripping section vapour flow, moles 


d(2,)qg = amount of any component i in the 
top product, moles 


b(z,;), = amount of any component i in the 
bottom product, moles 


Kiem and D. N. 


HANSON 


relative volatility of component 7. 
All relative volatilities to be based 
on the same component. 


¢ = a rectifying section parameter 
¢ = a stripping section parameter 


UNpDERWoOOD and (2) 


to define the rectifying and stripping section 
parameters, ¢ and ¢'. It is apparent that for 


employs equations (1) 


this three component system there are three 
values of ¢ which will satisfy eq. (1) such that if 
these are numbered in decreasing order of 
magnitude : 

a4 >>, > op 

ty > 2 > %& 


Xo > os 


There are also three values of ¢' which will 
satisfy equation (2) and if these are again listed 


in decreasing order of magnitude : 


o> 44 
a4 >o_g>%p 
tp > o's > && 

UnbeERWoop has shown that two independent 
equations may be written relating the three 
values of ¢, the concentrations of A, B and C 
in the liquid on the feed plate, and the number 
of rectifying plates. These equations are : 


Xp (Xp)y : a (Xe), 
xn—dy ac — dy by 
. Xp (2p), +4 ae (%e)y ¢; 


tp—o2 to-do 


ya) 


to (Le), 
to— de al , 
= 3 
te(te)y \de *) 
%co— ds 
where (x,), = mole fraction of component i in the 
liquid on the feed plate. 


X_ (Lp), 
. + 
an—de 
Xp (py 
Bes + 
an—dy 


n = rectifying plates, including the feed 
plate. 


The authors, to simplify the writing of such 
equations have adopted a shorter nomenclature, 
where 
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“A A,, “8 __. B,, ete. 
imme ¢, 2 | ?, 


and 
¢. 
$) 
Rewriting equations (3) and (3’), these 
become 


) = tes 


then 


A, (x4), + B, (rp), + 


Ag (@4)p + Bo (Xp)y 4 
and 


A, (24), 
As (@4)y 


C, (7%), 


, 3 
Cs (te) Yan 


+ By (@p)y + Cy (2), 
T Bz (@p)y + Cs (rc) 


two independent equations can be 


bs (3°) 


Similarly, 
written for the stripping section : 

A,’ (@ 4) oF B, (Tp)y + Cc,’ (%e) 

A,’ (T4) — B,'(®p)y T C;'(2c)s 

A, (@4)f — B,’ (7p), + C,’ (ty 

Ag’ (xv 4)s + B,’ (tp)y T C2’ (@e)y 


Hey 


: 32 (4’) 


where : 
% 
B sag 
2B dz 
$y )" 
>; 
where m = stripping plates, including the reboiler 
but not the feed plate. 
If components B and C are the key components 


‘te. 


and 


Yer 


and their separation is set, then the amounts 
of these components in the top and bottom 
products are known. As an approximation 
for the first solution, it will be assumed that 
component A goes completely into the top 
product. This approximation may be corrected 
to any desired degree of accuracy later and does 
not represent a loss in accuracy of the calculation. 
If, then, the reflux is set, the quantities V and 
V’ are known and the ¢ and ¢’ values may be 
calculated from equations (1) and (2). It will 
be noted that a consequence of the approximation 
of no component A in the bottom product is 
that ¢,°=a,4. This means that A,’ or 
a 4/(%4 — $,’) is equal to infinity and as a result 
eq. (4) becomes useless. Equations (3), (3’), 
and (4’) can be written out however and solved 
for n and m. 


Atper and Hanson have suggested that the 
solution of these equations be accomplished by 
solution of the third order determinant which 
can be set up as follows : 


Ag pay B, = B, Hey C, = 
Ay— Agta By— Bar Cy— Ca¥, 
A, - Ag $39 B, ¥ B, Psy Cs C, P32 


The and m are contained as 
exponents in the ~ and w’ terms. If a value of 
n is set, the corresponding value of m may 


C2 pa 


Cs bse 0 (5) 


unknowns n 


be calculated from eq. (5), and in this way the 
picture of all possible combinations of rectifying 
and stripping plates which will perform the 
distillation can be found. 

The solution of eq. (5) is not difficult. However, 
if, for example, a five component system were 
to be distilled, eq. (5) would become a fifth order 
determinant and the 
siderably more tedious. 


would be 
the number of 
components in the feed is equal to the order of 


solution con- 


Since 


the determinant to be solved, as the number of 
the 
However, by 


solution 
further 
development, the solution can be reduced from 
form to a 
form which is much easier to solve. 


components increases, labour of 


becomes prohibitive. 


determinantal difference product 


Kq. (5) may be written : 


| A, — Ag, Bid, — Bode Crh, — Code | 
Ay}, — Agts Baty — Baty Cat, — Cats | 
Ags — Ag}, Bybs — Bap. Cabs — C24. 


where 
py ¢,", Po = do", etc. 
bs = (1/b5)", Yo = (1/¢2)™, ete. 


It is apparent that the term a, a,%- can be 
cancelled out of the determinant, eq. (6). The 


ao 


term A, ¥,, for example, is then 


1 
a4— > 
Rewriting eq. (6) in a new nomenclature to show 
the dropping of all « values in the numerators 


of the individual terms is not worthwhile, 
however, and it can merely be understood that 
1 


for the rest of the development, 4, = ee 
ae 
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etc. This lapse in preciseness of nomenclature 
here leaves no ambiguity in the final solution. 

Eq. (6) can also be factored into six third 
order determinants to give: 


A, 4, A; 
B, B, B, 
C, C, C; 


A, A, 
Hy He vy B, By 
C.¢ 


na by bs 


A 
Vids, | B 
c 


If the first determinant in eq. (7) is written out, 
it is: 

1 
x4 bs | 





Such a determinant is a double alternant and 
can also be written [2] as 

(hk) (k#~1) 

(—1) 2% [8(a4 ap ae) f2(4, be bs) 


uy Up ic 


where & the order of the determinant. 


{? = difference product obtained by taking 
each term and subtracting from it 
in turn each term preceding, then 
multiplying all of these differences 

together. 

difference product obtained by multi- 

plying together all terms (a, — ¢,) 

where ¢, is a general value of ¢ which 
will take all values of ¢ used in the 
determinant. 

a difference product exactly like u, 

except that a, is used in place of «4. 

etc. 


Before proceeding, some simplification can 
be effected. Two factors enter into each deter- 
minant in eq. (7) and may be cancelled. Thus 
we may drop the terms 
(k) (k-—1) 
(—1) ° 
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and 


CF (a4 %p Xe) 
If these factors are dropped, then the first term 
in equation (7) becomes : 
(, te bs) ¢ (¢, be 3) 
U4 Upc 


and if this is written out it is : 


(yy de Wa) [($, — 2) (4, — 
[(x4 — $1) (a4 — 9) (a4 


$3) (2 — ¢s)) 
$s) 


[(ap — ,) (ap — oe) (a $5) 


[(xe — $1) (%¢ — $2) (ae — $5)] 

It will be noted that the liberty of writing the 
{*(¢, de 43) with all terms 
It is more natural to 


sub- 


difference product, 
reversed has been taken. 
calculate the difference product this way, 
tracting ¢ and ¢° values from ¢ or ¢ values 
preceding, and it yields the same results in the 
summation of terms in eq. (7). 

For illustration, the second term in eq. (7) 
Is 

(vy de He) 4 (by de oe) 


Uy Upc 
or written out 


(dy de He) [(F, 
[(a4 >) (%4 


bz) (¢, 
$2) (x4 


$2) ($2 — $,)] i 
$2) | 


¢;) (ap — $2) | 


1 
[(a¢ al ¢;) (a — $2) (a¢ — 


[(ap - $2) (%_ — 


$,)] 


The other terms in eq. (7) are constructed in the 
same way. 

At first acquaintance, the terms in the difference 
product solution of eq. (7) seem quite long and 
involved, and in the extension of the method 
to systems of several components, it appears 
at first that even the work of keeping track of 
the major terms in equations of the type of eq. (7) 
will be difficult. Since this is one of the chief 
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difficulties in the difference product solution, 
it is very important to set up a standard procedure 
and a set of rules that will make the difference 
product solution as automatic as_ possible. 
For this purpose, the general form of equation 
(7) may be written as: 


( 1) *Pte coe P+O+E+... 


UU, Ug. ..U, Us, Uy... 


(do Op q+ +r be Hees) 


Uy Uy Ug... U, Us, Uy... 


£4 (4, by by --- ob. $4) 
gcc Me Myeee | 

In general this equation is the solution of a 
determinant of the type of eq. (6). The 
order of the determinant is equal to the 
number of components in the feed, and in a 
case where the separation of the two key com- 
ponents has been set and the reflux has been set, 
one such determinant can be written [1]. This 
determinant will in general employ all of the 
¢@ values which are useful and all of the ¢’ values 
which are useful, thus excluding ¢ and ¢’ values 
which are equal to relative volatilities. The 
total number of useful ¢ and ¢’ values will equal 
the number of components plus two ; the number 
of ¢ values and the number of ¢’ values will vary 
depending on where the key components lie 





uu, u 


in the list of components. 
The rules which may then be set up for the 
solution of equation (7) are as follows : 


1. Each major term in the summation will be 
characterized by one less than the number 
of useful ¢ values and one less than the 
number of useful ¢° values. Thus the 
number of % values written down in each 


term is one less than the number of useful 
¢@ values, and the number of y values 
written down in each term is one less than 
the number of useful ¢’ values. 


The total number of ¢ and y’ values 
written down in each term is equal to the 
number of components. 


The y values are written down in ascending 
order of magnitude of subscripts (descending 


order of magnitude of ¢). The w’ values 
are written down in ascending order of 
magnitude of subscripts (descending order 
of magnitude of ¢’). 


All of the % values are written down first, 
followed by the y’ values, 


The combinations of ~ and w’ values 
used are all possible combinations under 
rules (1), (2), (8) and (4) which do not 
contain any repeating ¢ or ¢’ values. Thus, 
the number of combinations possible and 
also the number of major terms in the 
summation is equal to the number of 
useful ¢ values times the number of useful 
d values. 


The difference product [3 (¢, dy oq 

¢, ¢,¢,....) is obtained by subtracting 

each ¢ or ¢’ value from all ¢ or ¢’ values 

preceding it in the order of listing of the 

y and y’ values and multiplying these 

differences together. Thus {3 (¢, bp y+ +> 

«+ br be b+ ++) = [(bo — bp) (bo — be) - - - 

. (¢, ~ ¢,) (¢, 1 $,) ($ . ¢;) oe | 

[(, — 4) - «= (bp — $4) (bp — 4) 

(bp — $4) «+ =] [(bq — $1) (by — $) 
p-+-]--++[(¢, — $,)(¢, — $)----] 


The number of u values in each term is 
equal to the number of ¢ and ¢’ values 
used in the term. For convenience these 
have been grouped and characterized by a 
¢@ or ¢ value rather than by component. 
Thus u, = (a, — ¢:) (ap — ¢)) (x¢ — $,) 
(xp — $4) .» Ug = (a4 — $9) 

(ap — $2) (ac — $e) (ap — og)... ete. 


Each term in the summation is multiplied 
by —1 or +1 depending on the sum 
of the subscripts of all ¢ and ¢’ values used 
in the term. Since the choice is arbitrary, 
the convention has been adopted that if 
the sum of the subscripts is even, the term 
is multiplied by +1 and if the sum of 
the subscripts is odd, the term is multiplied 
by —1. 


The rules set down here completely define 





Gerearp Kien and D. N. Hanson 


the solution of the determinant and enable 
calculation of the column. Some study of them 
is necessary, however, before much facility in 
such solutions is gained, and to aid in this 
understanding of the method an exaniple of a 
distillation column separating a five component 
system has been solved for presentation here. 
Consider a five component system of A, B, 

C, D, E. The feed to the column is saturated 
liquid of the composition shown in the table 
below, and the relative volatilities based on the 
heavy key, D, are as shown : 

X»p a 

0-05 5 

0-25 2-5 

0-40 2-2 

0-20 1-0 

0-10 0-6 


Components C and D are the key components. 
Specifications on their separation are set so 
that 95% of C goes into the top product and 
90°, of D goes into the bottom product. It is 
assumed as a first approximation that the light 
diluents, A and B, go completely into the top 
product and that the heavy diluent, E, goes 
completely into the bottom product. Thus 
the products, based on one mole of feed, are : 


d(x)q (ap 
0-05 — 
0-25 _— 
0-38 0-02 
0-02 0-18 
— 0-10 


The minimum rectifying section vapour flow 
for this separation is calculated to be 1-206 
moles/mole of feed. An operating rectifying 
section vapour flow of 1-6 moles/mole of feed, 
or approximately 130°, of minimum, was chosen 
for the calculation. 

With these variables set, there is an infinite 
number of rectifying and stripping plate com- 
binations which will perform the distillation, 
and the answer desired is the combination of 
rectifying plates and stripping plates which will 
constitute a minimum of total plates, or will 
give the optimum feed plate location. Generally 


the minimum in the curve of total plates verses 
feed plate location is very flat, and a wide choice 
of feed plate location can be made with little 
effect on the total plate requirements. 

The first step in solution of the problem is the 
calculation of the ¢ and ¢’ values. This calculation 
is the solution of the following equations : 

For values of ¢ : 

1-6 = 
5 (0-05) 
5-—¢ 


For values of ¢’ : 


2-5 (0-25) 
25—¢ 


2-2(0-38) 1 (0-02) 


=a tea 7 


2.2(0-02) 1(0-18) | 0-6(0-10) 

o2-—-4 1-¢ ' 66—¢4 
Equations (8) and (9) may be solved by any 
of the iteration methods available for solution 
of cubic and higher order equations. Actually, 
trial-and-error solution is quite rapid and has 
the advantage that in the solution for ¢ or ¢’, 
the quantities 2, —¢,, etc., are calculated. 
Since these differences are 
convenient to have them tabulated as a result 
of this first calculation. 

Solving eq. (8), then, the following values of 
¢@ are found : 


— 26 (9) 


of use later it is 


It will be noted that a value of ¢, is listed. If 
the assumption is made that d(rg), = 0, then 
¢s = 4%, and may be listed as such although 
it is of no use in solving the column. 
Similarly, from eq. (9) the following values of 
¢’ are found : 
= § 
= 2-5 
2-2311 
= 1-118 
= 0-6284 


Here again values are listed for ¢, and ¢,. Under 
the assumptions of 6(z,), = 0 and b(r,), = 0, 
the values of ¢, and ¢, are respectively equal to 
a, and ag. These values of ¢, and ¢, are of no 
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use in solving the column, but are listed in order 
to clarify the system for subscript numbering 
of the ¢’ values. 

Three significant figures are being carried in 
the calculations at this point, although it is 
apparent that this necessitates having more than 
three significant figures in the ¢ and ¢’ values. 

The ¢ values to be used in the solution are 
¢;, de, $3, and ¢,. The ¢’ values to be used are 
ds, $4, and ¢,. Generally, the next step in calcula- 
tion is to form a table listing the characterizing 
y and ¥’ combinations. In the table may be 
listed, then, the other terms which accompany 
the characterizing ¥ and w’ combination in 
equation (7). These are: in the numerator, 
the term (— 1) @*? 74+ --- 7 t# +" *-9. which will 
be called simply the sign; in the numerator 
the value of the difference product, {3 (¢ ); 
and in the denominator the term u,u,... 
Such a table has been constructed for the 
example and is shown in Table 1. In the first 
column are listed all of the characterizing ¥ and w’ 
combinations which are possible under the rules. 
As an illustration, the various items in the first 
row of the table will be calculated out here : 


The sign by which the term is to be mul- 
tiplied is determined by (— 1) ‘* ‘scripts on ¥ 
and # values) which, for the first combination, 
is (— 1) +#+8+*+9— (_ 1)'5— —1, The 


sign is thus minus. 


Next the difference product is, when 
written out : 


ti(¢ se) = [(¢, = $2) (¢, rz $s) (¢, rs $s) 
(d, ‘sar $s) [($z = $s) ($, 4 $s) 
(¢2 — $s) | [(¢s a $s) (d5 i ¢s) | 
((¢. — ¢s)] 
= [(2-49) (8-50) (3-77) (4-26)] 
[(1-01) (1-28) (1-77)] [(-264) (-754)] 
[(-490)] 
= + 81-0 
Lastly the product of u values = wu, uy uy 
u,u, where 
u, = (a4 — ¢) (4p — 4) (a — 4) 
(ap — $,) (ag — ¢,) 
= 12-3 etc. 


Then u, u, us uu, = + -0787. 


The various values of « — ¢, are available 
from the calculation of ¢, and it is convenient 
following calculation of the ¢ and ¢’ values 
to calculate all values of u and tabulate 
them for use. 


These terms are entered in the table, and the 
rest of the table is calculated, whereupon eq. 
(7) may be written : 


0-0787 


— $1 te ds Ha os (or) 


+o, be bs oa vs (Fis) 
— 4.54 
0-0269 

— 13-2 
ns ~onaé) = 


— oa da dn tn He ( 


+ th ba a da ds | 


If, now, a value of n, the number of rectifying 
plates, is set, from eq. (7) the corresponding 
value of m may be calculated. In this example 
the following values have been calculated : 


n-+-m 
10-5 
10-2 
10-5 


n m 
6 45 
5 5-2 
45 6-0 


It is apparent that the minimum in the curve 
of total plates is flat, and that n = 5 and m = 5-2 
essentially represents the optimum point. 

The assumption, made at the begining of 
the example, that components A and B, which 
are lighter than the light key, go completely 
into the top product, and that component &£, 
heavier than the heavy key, goes completely 
into the bottom product, should be checked. 
The actual amounts of these components in the 
products will vary to some extent with the number 
of rectifying and stripping plates. Since the 
minimum point in total plates is of most interest, 
the amounts of the diluents in the products 
should be calculated at this point, i.e., with 
rectifying plates equal to 5 and stripping plates 
equal to 5-2. The calculation of the amount of 
each diluent in the products requires the solution 
of one additional determinant. This is identical 
to the solution for plates, except that the ¢ 
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and ¢° values used are different. The correct 
values of ¢ and ¢’ to be used for the checking 
of each diluent is discussed in the previous 
paper by ALper and Hanson [1]. Following 
the calculation of the amounts of the diluents 
in the products, the question of whether or not 
the calculation for the number of plates should be 
repeated with more correct ¢ and ¢’ values is a 
matter of judgment and depends on the accuracy 
desired in the final answer. 

The calculation of this example, while it may 
appear so at first glance, was not excessively 
time consuming. Actually the most labour in 
such a calculation occurs in calculation of the 
¢@ and ¢’ values. Familiarity with the rest of 
the calculations leads to many short cuts which 
save considerable time. For example, it should 
be apparent that if eq. (7) is multiplied through 
by the product of all u values, there remain 
in each major term only two u values in the 
numerator, the values of u corresponding to the 
one ¢ value and the one ¢’ value which do not 
appear in the characterizing y term. Along with 
this simplification, others can be effected, and 
convenient tables can be set up for calculating 
the difference product terms. The amount of 
time which will be spent on the calculation of a 
column by this method is a definite function of 
the amount of experience with the method which 


the calculator possesses. Whether the method 


Table I. 
Characterizing | 
term «Sign fi(¢....) | 


Sample calculation 





Uy Up... 





+ 310 + 0-0787 
— 29-2 — 0-0103 
4°54 + 0-0269 
— 13-2 0-00956 
30-4 0-00125 
6-64 0-00328 
0-484 0-0705 
28-4 0-00922 
3-00 — 0-0241 
0-00717 0-0007 67 
0-0758 0-000101 
0-00655 0-000263 


tr do Hs He os 
+; ¥2 43 3 Hs 
+, ve vs Hy He 
Hy be ve He Hs 
+, ve vs bs vs 
¢, ete da ve + 
Ys Vededs | - 
d, os He Us Hs 

ty os Ye Ys Me 

te os 5 Ue Us 

te oy He ¥y Vs 

be ty be Hy Me 





is efficient in comparison with plate-by-plate 
procedure is mostly dependent on the system 
to be fractionated, and again it must be left 
to the judgment of the calculator to determine 
whether it is desirable to employ this method. 
Certainly there are many cases where it will 
prove of considerable value. 


NOTATION 


amount of component i in the bottom product, 
moles /mole of feed. 
= amount of component i in the top product, 
moles /mole of feed. 
order of a determinant. 
rectifying plates, including the feed plate. 
stripping plates, including the reboiler. 
product of all terms (x4 — ¢,), where ¢, is 
any value of ¢ or ¢' occurring in the alternant 
being solved. 
= similar product of all terms (ag — ¢,), ete. 
product of all terms (a; — ¢,), where a; is 
the relative volatility of a general component. 
The product contains a term for each com- 
ponent. 
product of all terms (a; — $y). ete. 
mole fraction of component i in the liquid 
on the feed plate. 
etc., refer to components ; in general A being 
the lightest, B next lightest, etc. 
4, B, —_ 
x4— tp? 
‘ , , ae 
= abbreviation for ; B,=- 33 
24-9) a 
* = rectifying section vapour flow, moles /mole of 
feed. 
)’’ = stripping section vapour flow, moles/mole of 
feed. 
(X,)p = mole fraction of component i in the tota 
feed. 


= abbreviation for etc. 


etc. 


Greek Letters 
a, = relative volatility of component i. All 
relative volatilities are based on the same 
component. 
rectifying section parameter. 
stripping section parameter. 


¥%, = abbreviation for ¢,", ¢e = 2", ete, 
¢, = abbreviation for 1/(¢,)", vg = 1/(¢g)", ete. 
; : b, n 5 n 
¢g, = abbreviation for [—"] . ¢g2 = ( , ete. 
¢) +, 


a \m 
_ ) , ete. 
%, 


"\m 
$y, = abbreviation for (*) » 32 = 
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{i = difference product obtained by subtracting General Subscripts 


each term in an alternant from each preceding 0, p,q... 1’, 8’, t’ refer to numeral subscripts on ¢ and ¢’ 
term and multiplying these differences values. 


together. 


i refers to a general component. 
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The viscosity temperature relationship of liquids 
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Summary—aA new formula representing the viscosity temperature dependence is proposed. 
The formula is valid for a variety of liquids, e.g. associated and non-associated liquids, hydro- 
carbons and mixtures of hydrocarbons. Water, solutions of sucrose in water and of sulphuric 
acid in water also fit into the system. The same formula can even be applied to mercury. The 
method can be used for saturated and non-saturated natural mineral-oils, for fatty oils as 
well as for hydrocarbon mixtures obtained by polymerization and by expelling oxygen from 
polymerized fatty oils. The proposed formula is generally applicable and more concise and 
clear than many of the existing equations. 


Résumé—Pour représenter la relation température — viscosité, les auteurs proposent une 
nouvelle formule valable pour une série de liquides : liquides associés et non associés, hydro- 
carbures, mélanges d’hydrocarbures, eau, solutions de sucre dans l'eau et d’acide sulfurique 
dans l'eau, mercure, huiles minérales saturées et non saturées, huiles grasses, mélanges d’hydro- 
carbures obtenus par polymérisation et par réduction d’huiles grasses polymérisées. La formule 
est générale, plus concise et plus claire que beaucoup d’équations existantes. 


INTRODUCTION % 


80 60 





It was the intention of the authors to contribute 
to the study of the connection between viscosity 
and change of viscosity with temperature res- 
pectively and the structure of hydrocarbons and 
mixtures of hydrocarbons. The first results in 
this field were obtained with liquid polymers 
(from definite monomers, olefins of more or less 
branched character), sometimes containing 
rings [1, la]. Regularities were noticed, parti- 
cularly in the temperature dependence of the 
viscosity. Use was made of the de Guzman- 
Andrade formula [2, 3], in this case applied to s 
mixtures of substances : ry 






































A $00) 
logy =< +B ie 


7 = dynamic viscosity in cp. 











l = 25 ise : j _ L aemeee 
P= ars pe = A 
T = absolute temperature °K. 


A and B = constants. 





























This formula could be used for classification : 7300 3500 
purposes, but deviations occurred with fractions co on 
of higher molecular weight (Fig. 1) beginning with : 
fraction P-83,, viz. to begin with the fraction Fis. 1. log» — 5 diagram for polymer-fractions P-1, 
having a value of 747°°© equal to about 1-5 cp. P-2, ete. 
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The investigations reached this point before the 
war and were resumed after the war. New results 
have been obtained, which have partly been 
published [4, 5, 6]. 

It must be remarked that the study of the 
viscosity temperature relationship is not new : 
other investigators tried to find formulas con- 
necting viscosity with temperature. 

One of the first attempts has been made by 
PoIsEUILLE, who gave the formula [7] : 


7 = (1 + of + Bt) (2) 
» = dynamic viscosity in cp. 
% = dynamic viscosity at t = 0°C, 
t = temperature °C, 
aand § = constants. 


Much use is made of the de Guzman-Andrade 





















































formula (1); see also [8]. Soupers [9] replaced 
dynamic viscosity by kinematic viscosity and 
obtained : 


logy = +B (3) 


1 
= kinematic viscosity in cs. | 1 cs = —— 
v ¢ viscosity in c ( cs 100 


2 
stokes ; stokes = =). 
sec 


T = absolute temperature K. 
A and B = constants. 


Both formulas (1) and (8) fail with associated 
liquids and with the hydrocarbon mixtures of 
higher viscosity discussed above, as was already 
known in the literature for lubricating oils [14]. 

WaLTHER [10] proposed the formula : 


log log (v + a) = W — mlog T (4) 


c 
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1 
Fic. 2. log » — 7, diagram for a number of hydrocarbons. Fie. 8. log » — 34 diagram for s number of hydrocarbons. 
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v = kinematic viscosity in cs. 
T = temperature °K. 
a, m, and W = constants. 


This formula is often applied to lubricating 
oils. For general use the difficulty is to determine 
the value of the constant a. For lubricating oils 
the value of a is often assumed to be approximately 
constant, and was considered to have an average 
value which has been reported successively as 
0-8 and 0-6, though a better approximation 
appears to be 0-7 [10a]. 

Voce. [11] proposed the formula : 


(5) 


™ = dynamic viscosity in ep at °C. 
1. = dynamic viscosity in cp at t = 0. 
t = temperature in °C. 
t, = temperature at which », = 1. 


t_. = temperature at which n, = ©. 


The formula may be written as : 


B 
logn=A+FTe 


» = dynamic viscosity in cp. 
T = absolute temperature °K. 
A, Band C = constants. 


In this form the formula is discussed by Gut- 
MANN and Simmons [12]; see also [13, 14]. 
°C 
9 8 1 60 
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Fic. 4. logy — 7s diagram for a number of hydrocarbons. 


Recently another formula giving the viscosity 
temperature dependence was proposed by 
UmstAtrer [15], which may be written as: 


are sh In = B+ Cin T (7) 
7 = dynamic viscosity in cp. 
T = absolute temperature °K. 

A, Band C = constants. 


This formula is discussed by Frssier [16]. 


THe New FormMuta 


It appeared to the authors that the deviations 
from the de Guzman-Andrade equation (1) or 
rather the SoupErs equation (8) were very 
regular (see Fig. 1). 

A new formula, resembling that of Soupers, 
which is hereby proposed, appears to hold in 
cases in which the SoupERs equation was not 
applicable : 


(8) 


A 
l =— +B 
og v jet 









































Fic. 5. logy — 7a dingram for water, H,SO,-H,O 
mixtures and sucrose solutions in water. 
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v = kinematic viscosity in cs. 
T = absolute temperature °K. 
A, Band z = constants. 


Graphically this equation is represented by a 
straight line if log v is plotted against a for a 


T 
particular liquid; B is the value of log v at 
T= 0; A gives the temperature dependence 
of the viscosity of the liquid under examination. 

The formula was applied to many liquids ; 
single liquid compounds as well as mixtures of 
single liquid compounds, non-associated as well 
as associated liquids. 

In order to calculate the value of z in equation 
(8) for a liquid we have to make three viscosity 
measurements, v,, vg and vg at three tempéra- 
tures T,, T, and 7;, respectively, 

°C 


20 40 
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- 
Fis. 6. log » — 74 diagram for hydrocarbon mixtures 


obtained by removing oxygen from polymerized fatty oil 
fractions. 


A 
log = pe t+ B 
2 


A 
From the former two equations it follows by 
subtraction that 
er 
6 "1 — 06 "9 = Te: Te (a) 
and from the latter two : 
log vg — log vs = we 


(a) and (b) are two equations with the two 
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Fic. 7. logy — 7a; diagram for a number of lubricating 
oils. 
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unknown constants A and z. By division, the 


equations can be transformed into : 


_ logy, — log ry | T 
~ log vy — log vs 


From this equation z may be calculated if 
Vy, Yq. Vg» Ty, T, and T, are known. This can be 
done by calculating the value of Q for various 
values of z. 

As the viscosities were determined at 20°C, 
40°C and 70°C (T, = 293°K, T, = 313°K and 
T, = 343°K), the value of Q may be written as 


20 
Q< 40. 
70 


The results are reported in Table 1. 


On the other hand if the viscosities v,, v, and vs 
are determined at the temperatures 20°C, 40°C 
and 70°C, the value of Q may be calculated from 
equation (9a) : 


log Vy 


log vg 


log vg 


— log vs 


Q (9a) 


By means of Table 1 the value of 2 may be 


of Q. 


determined from the calculated value 


Table 1. 





20 

iF) wy 

70 
0-78 
0-84 
0-91 
0-98 
1-06 





An important feature of the new formula proves 
to be that the nearest integral value may be 
assigned toz without introducing excessive errors in 
graphical interpolation or extrapolation provided 
the temperature range is not too large. E.g. the 
determined value of z= 2-2; it now proves 
permissible in equation (8) to use z = 2-0, instead 


of z = 2-2. This gives only a slight deviation 
from the straight line in the diagram log v = f 


1 
(7) in which 2’ = nearby whole value of z. It 


means that only 4 or 5 diagrams (for z = 1 to 
a = 5) are necessary for practical use. It also 
means that for the various members of a group 
of liquids the values of A, the constant which 
indicates the temperature dependence of the 
viscosity, may be compared with each other with 
the aid of the nearest integral value of 2. 


EXAMPLES 
A large number of viscosity values were deter- 
mined by means of the Wi. OstwaLp viscometer 
with improvements worked out by Sproxet [17]. 
Other viscosity values were taken from literature 


c 
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Fic. 8. logy — Th diagram for a number of lubricating 

oils. 
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The viscosity temperature relationships of liquid 


Table 2. 


| | 


Hog "20 — 10R ¥49 _ 
Compound log v¥49 — log v79 (¥30)pred. (¥g0)exp. A% | (¥40)pred. (¥50)pred. (so)exp. 











0-886 | 0-885 0-1) 0-795 -4| 0-719 | 0-716 
0-895 | 0-885 11| 0-794 8| O712 | 0-716 
n-Pentadecane . 8-05 0-7; 2-523 -2| 2-132 | 2-128 
n-Hexadecane . 2-954 : 2-473 | 2-479 


n-Nonane 








Et-Cyclopentane . 0-609 0-608 ° 0-558 | 0-557 
n-Pentylcyclopentane . : . 1-103 1-102 | -1| 0-978 | 0-978 





Cyclohexane | 0-924 | 0-926 2) 0-808 | 0-808 
n-Butylcyclohexane . | 1-409 . . 1-222 1-217 : 1-074 | 1-069 
n-Tridecyleyclohexane ;, 8-18 9| 627 | 6-26 2| 492 | 4-99 
n-Decyleyclohexane , 5-01 6, 397 | 3-96 8) 321 | 3-24 


n-Nonene-1 . 0-762 | Oo 5| 0-688 | 0-690 |—0-3! 0-624 | 0-625 
n-Pentadecene . 2-64 , —O-4 | 221 2-22 5) 1-89 | 1-89 








n-Propylbenzene . 0-8741 | 0-8745 | 07759 | 0-7770| — 0-1) 0-6972 | 0-6979 
n-Butylbenzene . | 1048 | 1-049 1) 0-923 | 0-925 | ‘2| 0-822 | 0-818 
n-Hexylbenzene 1-661 | 1650 | 08) 1-426 | 1-419 ‘5| 1-240 | 1-239 
n-Heptylbenzene . 2-04 2-03 ‘S| 1-722 | 1-722 0} 1-483 | 1-488 
n-Undecylbenzene . | 428 423 2); 345 | 343 | 6| 2-84 | 2:85 
n-Pentadecylbenzene . 780 | 792 |— 15) 605 | 615 6) 483 4-92 

















(s0exp. | A% 





- | 0-552 “ | © “ 0-510 
n-Nonane 6 | . ’ 0-542 . . -f 0-510 


n-Pentadecane 2) | . 1-393 : 5 | 12% 1-244 
n-Hexadecane " . . 1-583 , 6) 1 1-408 


Ei-Cyclopentane , . | 2| 0442 | 0 ‘5 . | O41 
n-Pentyleyclopentane ’ | 0-722 : : . | 0-67 


Cyclohexane : . 2 0-569 
n-Butylcyclohexane , . . ‘6 | 0-776 
n-Tridecylcyclohexane , 3) 2-81 
n-Decylcyclohexane . . , . 1-95 


n-Nonene-1 
n-Pentadecene 


n-Propylbenzene — 01 
n-Butylbenzene 0-5 
n-Hexylbenzene 0-1 . , . 0-876 
n-Heptylbenzene — 0-4 , ‘7 | 1024 
n-Undecylbenzene — 0-4 ‘5 | 1-80 

n-Pentadecylbenzene | —18 : 3-36 2-83 | 
































In this table » is given in centistokes. 
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data. For each liquid the value of z was calcu- 
lated by means of equation (9a) and Table 1. 

In the first place the results of interpolation 
of viscosities at various temperatures of some 
hydrocarbons [18] will be shown. The values of 


20 

Q (< « and z are calculated from the 
70 

viscosities at 20°C, 40°C and 70°C. The viscosities 

at 20°C, 60°C and 100°C were used to construct 


1 
the straight line in the diagram log vy = f (7): 


The viscosities at 30°C, 40°C, 50°C, 70°C, 80°C 
and 90°C were predicted and compared with the 
experimental values from literature, see Table 2 
and the corresponding diagrams, Figs. 2, 3 and 4. 
Table 2 shows that the nearest whole value of 
aw may be taken without making large errors in 
the interpolation. E.g. the value of z of n-nonane 


‘| 
is 1-5 (Table 1). Either the diagram log vy = f (7) 
1 
or the diagram log v =5(p) may be used. 
The consequent deviation rarely exceeds 1-5%, 
which is about equal to the experimental error. 


The nearest whole value of z may be taken only 
if the temperature range is not too large, e.g. 
from 20°C to 100°C. If the range covers several 
hundreds of degrees Centigrade a more exact value 
of z has to be chosen (see below). 

Further examples are given in Table 3. See 
also Fig. 5. 

In Table 3 predicted values (obtained by 
graphical interpolation) and experimental values 
are given of a sucrose solution in water (con- 
taining 60% sucrose by weight), a H,SO,-H,O 
mixture (containing 20% H,SO, by weight) and 
water itself. 

20 

The value of Q (< { | of the H,SO,-H,O 
60 

mixture is calculated from viscosities at 20°C, 

40°C and 60°C. The corresponding value of z 
can be read from Table 4. 

As can be seen from Table 3 good agreement 
exists between predicted and experimental values 
of viscosity. 


In Fig. 5 data for water, two H,SO,-H,O 
mixtures and three sucrose solutions in water are 


Table 8. 





60 %, sucrose-40 %, H,O 


20 % H,SO,-80 % H,O 





(Yes pred. 4% 


(%cg exp. 


(ma exp. 4% 





SSSSSESESSS 


a 











23-5 0-9 
17-55 0-2 
13-56 0-7 
10-80 12 

8-04 0-1 


6-48 0-0 





























The viscosity temperature relationship of liquids 








1-29 


1-46 








1 
plotted in a log vy — T° diagram, T ranging from 


273°K-873°K (0°C-100°C). 
In Fig. 6 data for mixtures of hydrocarbons 


obtained by removing oxygen from polymerized 


' 1 
fatty oil fractions [19] are plotted in a log v — i 


diagram. Measurements were made at 20°C, 
40°C and 70°C. The value of z for all these 
fractions appeared to be about 3. The viscosity 
values ranged from 5cs to 700cs. The hydro- 
carbon mixtures contained from about 0 to 3 
rings per molecule. 

Next data of 25 lubricating oils are plotted in 


1 
two logy — a diagrams (Fig. 7 and Fig. 8). 


Measurements were made at 20°C, 40°C and 70°C. 
20 

40 for these lubricating oils 
70 

ranged from 0-91 to 0-99 and consequently the 
values of z from 8 to 4 (see Table 1). The average 
value of z = 34 was chosen to construct the 


The values of Q 


1 
log v — rs diagram. 
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1 
Fic. 9. logy — Tih diagram for mercury. 
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Finally data of viscosity values of mercury [20] 
are shown in Fig. 9. The value of z appeared to 
be z = 1}. The temperature ranges from 0°C to The authors intend to develop their research 
310°C. As may be seen from the diagram work in various directions. 


1 
logy = f (yn) mercury also satisfies equation (8). 
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Book reviews 


Chemical Engineering in Practice. Edited by James 
I. Harrer. Reinhold Publishing Corporation. 24s. 


In May, 1953, a meeting was held by the American 
Institute of Chemical Engineers at the University of 
Pennsylvania at which eight authors read papers con- 
cerned with the putting of a chemical process into 
successful commercial operation. These papers have been 
edited by J. I. Harper and are collected in this small 
volume to put before the young chemical engineer some 
of the problems that may beset him and the way they 
can be solved. 

The first three chapters are on Process, Research and 
Development engineering and are illustrated by the 
laboratory investigations in the production of high 
octane gasoline from catalytic reforming of petroleum 
naphthas ; the design and construction of pilot plants, 
and the selection of the most economic design illustrated 
by the production of hydrogen, either by hydrolysis or 
by propane reforming. 

The next chapter is on Economic Analysis and shows 
how the annual return on investment or the payout 
period (i.e. the number of years that will elapse before 
There 
are then three chapters on Project, Construction and 


the invested cash is recovered) can be estimated. 


Operational Engineering corresponding to the problems 
encountered in detailed design, in construction on the 
site in operating the plant to bring it into production. 

The final chapter is on * Market Research Engineering ” 
although it is difficult to see how the term engineering 
ean be included. It is illustrated by those amateurish 
pseudo-comic sketches beloved by the Americans, but 
which only raise a puzzled expression on the faces of their 
foreign colleagues. Market research has so many impon- 
derables that it would be a rash manufacturer who 
entrusted this work to junior and unexperienced staff 
to whom these lectures are addressed. 

It would have been more instructive if the eight authors 
could have collaborated to show these stages of develop- 
ment with a single process as was attempted in “A 
Problem in Chemical Engineering Design” issued by 
the Institution of Chemicai Engineers for the more 
limited field of process and project engineering. The 
price is very high compared with other books now available 
in the chemical engineering field. 

M. B. Donato 


Warren K. Lewis, Artuur H. Rapascn, H. Ciay 
Lewis : Industrial Stoichiometry. McGraw-Hill Book 
Co. Inc. 35s. 


Tue sub-title of this book is “ Chemical Calculations 
of Manufacturing Processes.” It was first published in 
1926. The dictionary definition of stoichiometry is “ the 


branch of chemistry treating of chemical combination 
in definite proportions.”’ To those who took the chemical 
engineering practice school at M.I.T. it was the method 
by which the efficiency and method of operation of any 
chemical process could be studied without the necessity 
of the paraphenalia of weighing machines, flow meters, 
tapping holes in reaction vessels, the means by which 
the process could be analysed with the least possible 
interference with its normal operation. It is to be hoped 
that one day mechanical engineers will realize its usefulness 
in assessing boiler performance and drop the calculations of 
weight of air required for combustion when they invariably 
order their fans on volume. 

The change in the second edition is indicated by the 
increase in size from 174 to 429 pages. The chapter on 
Fuels and Combustion is split into Gaseous and Liquid 
Fuels on the one hand and Solid Fuels on the other. The 
chapter on Gas Producers is now headed Secondary Fuels. 
There are additional chapters on Energy Balances and 
Equilibrium, Nitrogen Compounds, Fixed Alkalis and 
Ceramics. Crystallization and furnace design chapters 
disappear. The chapter on Equilibrium contains a com- 
prehensive graph of the equilibrium constants for com- 
bustion reactions against temperature; Nitrogen 
Compounds include ammonia and nitric acid, Fixed 
Alkalis deals with the manufacture of soda ash and its 
causticisation ; Ceramics treats the use of binary phase 
rule systems and that of magnesia-aluminia-silica for the 
solution of problems. 

If the book is not already the standard text book of 
all schools of chemical engineering it certainly should be, 
as indeed it can almost now be included under the category 
of a classic in the subject. 

M. B. DonaLp 


W. Martz (Farbwerke Hoechst): Anwendung des 
Aehnlichkeitsgrundsatzes in der Verfahrenstechnik . 
Springer-Verlag, Berlin, 1954. 115 pages. DM 13.50. 


THe main purpose of this book, as stated in the preface, 
is to introduce the reader to the principle of similitude and 
to the theory of models, as well as to show the possibilities 
of practical application. After an introductory chapter 
entitled “General Remarks’ the author in the second 
chapter deals with the “ Theory of Similitude.” It is a 
difficult task to explain in an exact, logical and at the same 
time clear and simple form this theory and the theory of 
dimensions. far more difficult than it is to apply the results 
of these theories to practical problems. Dr. Matz does not 
intend to contribute to the theoretical aspects. It is a pity, 
therefore, that in trying to be complete, he develops the 
theory from three different starting points, viz. from the 
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proportionality factors of the physical constants, from the 
differential equations governing the phenomena and from 
the physical dimensions of the quantities under study. For 
the reader who studies this subject for the first time this 
multitude of derivations may be embarrassing. Besides, 
the treatments are not quite exact. E.g., the foundatien 
of the theory of similitude on differential equations re- 
mains rather vague. It does not become clear from the 
treatment that the differential equations of two similar 
systems become identical if both the coefficients and the 
boundary conditions are made dimensionless. The author 
has not made use of the opportunity to elucidate in this 
section the definition and the implications of the principle 
of similitude of two systems from a mathematical point 
of view. 


The third and most extensive chapter of the book is 
devoted to engineering applications. In the examples the 
author gives rather a large part of the normal non-dimen- 
sional theory of the unit operations involved and from 
this the reader is led to the rules for similarity and the use 
of dimensionless numbers. In a book that is intended 
to explain the methods of dimensional analysis it would 
have been more logical to stress the fact that these methods 
can be used without solving the problem analytically. 


In some of the examples, where the author shows that 
exact similarity is not possible, an explanation is given of 
very useful artificial corrections which can be applied to 
the model in order to reach a satisfactory degree of 
similarity. 

Some details of the treatment may be criticized. The 
reviewer cannot agree, for instance, with the treatment 
of the power consumption of a paddle stirrer on page 57, 
where the power number is given as a function of the 
Froude number only. Moreover, the power number as 
plotted on the ordinate of Fig. 12 is not dimensionless. It 
is the reviewer's personal feeling that some other intro- 
ductions to the subject which appeared in the English 
language (e.g. the excellent book of L. Lancuaar, Dimen- 
sional Analysis and the Theory of Models) are to be pre- 
ferred, but for those to whom these English texts are not 
accessible this book provides a sufficient basis for the 
application of the very useful principles of similitude. 

Th. N. Zwrererino 


H. Fiscner: Elektrolytische Abscheidung und 
Elektrokrystallisation von Metallen, 717 pp. 247 illus- 
trations. Springer, Berlin. DM. 72.00 


In us preface to the book Professor Fiscuer states that 
it was his purpose to write a fundamental (grundlegendes) 


work covering the whole range of scientific principles con- 
nected with the study of the electro-deposition of metals. 
In so doing he hoped to establish a scientific foundation 
for the relatively new subject of “ electro-crystallisation.” 


The first part (some 320 pages), of the book is devoted 
to the fundamentals of electrochemistry with particular 
reference to cathodic processes. It includes a section of 
over one hundred pages dealing with the role of inhibitors. 


Part two deals with “ electro-crystallisation.” After 
introductory sections on the kinetics of crystal growth in 
general and on the mechanism of electrolytic crystallisation 
in particular, the author goes on to discuss the structure 
of electrodeposits, presenting his material in terms of the 
system originated by him some five years ago. There is 
a short final section dealing with non-metallic inclusions. 


Part three discusses the theoretical aspects of the most 
important properties of electrodeposits, including bright- 
ness, internal stress and hardness. The forty-two pages of 
part four list values of some of the above properties for 
the deposits of many of the important metals. The final 
five pages are devoted to alloy deposits. 


The reader of this book will find that in it is collected 
a large range of material which has hitherto been dispersed 
throughout the literature. The book is well indexed and 
it will therefore prove valuable as a source book, although 
the text itself is not designed for easy reference. 


It is unfortunate that a tendency to over-elaboration 
makes the book difficult to read. All the data 
relating to the structure of polycrystalline electrodeposits 
are arranged in terms of the previously mentioned system 
of the author’s. One wonders whether this part of the 
book would not have been more useful as a source of 
information, had the facts been presented independently 
of the system and the latter described in a separate 
section; it would have benefited by more thorough 
illustration and the microphotographs could perhaps have 
been selected from a wider field. One would, for example, 
have expected to see one of the classical photographs of 
Erpey-Gruz and Vo-mer of silver or copper deposits on 
single crystal spheres. Many readers will be a little 
worried by explanations in terms of, for example, “ streng- 
thened, but still moderate, inhibition.” 


These criticisms having been made it is necessary to 
state that Professor Fiscner has rendered an invaluable 
service by undertaking such a laborious task of compila- 
tion. There is little doubt that the book will prove valuable 
to anyone concerned with the fundamentals of electro- 
deposition. 

H. J. Pick 





